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1.0  PASSIVE  SONAR  PERFORMANCE  ANALYSIS 

By  the  end  of  World  War  II,  available  methods  for  predict- 
ing and  analyzing  the  detection  performance  of  sonar  systems 
were  well  documented  in  the  Underwater  Sound  series  of  the 
summary  reports  produced  for  the  National  Defense  Research 
Committee.  These  methods  can  be  summarized  by  two  entities: 
the  sonar  equation  and  the  receiver  transition  curve. 

The  objective  of  this  report  is  to  develop  methods  for 
predicting  the  performance  of  passive  sonars  operating  in 
environments  with  fluctuating  sonar  equation  parameters. 
Although  this  subject  has  received  considerable  attention 
in  the  last  two  decades,  nearly  all  of  the  effort  pertains 
to  cases  in  which  the  relaxation  time  of  the  fluctuations  is 
large  compared  to  the  post-rectification  averaging  time.  The 
end  product  of  the  development  is  a set  of  transition  curves 
applicable  to  a limited  set  of  fluctuating  acoustic  environ- 
ments. The  use  of  these  curves  requires  evaluation  of  this 
sonar  equation  to  determine  the  average  excess  signal-to-noise 
level  for  particular  sonar  operating  environements. 

The  remainder  of  this  section  discusses  the  rationale 
and  utility  of  the  sonar  equation  and  the  transition  curve, 
and  identifies  the  essence  of  the  current  state-of-the-art. 

This  discussion  serves  as  an  introduction  for  a description 
of  certain  salient  features  of  the  approach  followed  in  the 
remainder  of  the  report,  and  also  serves  as  a framework  for 
the  utilization  of  the  resulting  receiver  transition  curves. 
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A block  diagram  of  a broadband  multibeam  sonar  system  is 
given  in  Figure  1.1,  in  which  the  heavy  arrows  indicate  multi- 
ple paths.  Only  one  of  numerous  sources  of  background  noise 
is  represented.  The  prediction  problem  can  be  divided  into 
two  parts,  as  suggested  by  the  vertical  dashed  line.  To  the 
left  of  that  line,  there  are  components  that  are  generally 
regarded  as  having  linear  transfer  characteristics,  such  as 
the  acoustic  transmission  channel,  the  sensor  elements  of  the 
array,  and  the  bandpass  filters  of  the  receivers.  The  beam- 
forming operation  may  be  linear  or  nearly  so,  or  be  non-linear 
when  the  inputs  are  hard-clipped.  If  all  of  the  components 
are  linear,  then  the  concept  of  signal  power  and  noise  power 
at  the  outputs  of  the  bandpass  filters  is  unambigious.  The 
principal  difficulty  here  is  predicting  the  characteristics 
of  the  acoustic  transmission  channel  from  the  source  of  interest 
to  the  array. 

A variety  of  processing  functions  are  performed  in  the 
portion  of  the  system  on  the  right  side  of  the  dashed  line. 

Each  receiver  includes  a rectifier,  a device  with  an  inherently 
non-linear  transfer  characteristic.  The  outputs  of  the  recti- 
fiers are  shown  on  an  intensity  — or  reflectance  modulated 
display  in  a beam-time  format.  An  operator  views  the  display 
from  time  to  time  and  makes  decisions  regarding  the  presence 
of  the  source  of  interest.  The  analysis  of  operator  performance 
is  considered  a very  difficult  undertaking;  however,  recent 
researches  in  neurophysiology*  have  produced  information  regard- 
ing important  functions  performed  in  the  human  visual  perception 
system. 
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•Reviewed  in  Section  2.0  of  Reference  4. 
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FIGURE  1.1  Representation  of  Multibeam  Sonar  System. 
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Given  these  complexities,  it  has  been  found  expedient  to 
consider  the  elements  to  the  right  of  the  dashed  line  as  a 
unit,  and  to  characterize  their  combined  performance  by  means 
of  a transition  curve  that  can  be  obtained  experimentally.  A 
transition  curve  gives  the  probability  Pg  that  the  operator 
will  detect  the  presence  of  the  source  of  interest  with  one 
observation  of  the  display,  as  a function  of  the  difference 
Ng  of  the  average  signal  level  and  the  average  noise  level 
that  would  be  measured  at  the  outptut  of  a bandpass  filter  of 
the  beam  directed  at  the  source  of  interest.  The  value  Np  of 
Ng  for  which  Pp  equals  a specified  value  is  called  the  detec- 
tion differential,  or  the  recognition  differential  if  classi- 
fication of  the  source  is  achieved.  Usually  the  value  specified 
for  Pp  is  0.5.  The  assymptotic  value  for  Ng  -►  -00  is  Pp,  the 
probability  of  false  alarm  for  a single  observation.  If  the 
operator  has  an  option  regarding  the  duration  of  record  displayed,  • 
then  the  chosen  option  should  be  specified. 

A transition  curve  can  be  plotted  as  a function  of  Ng,  or 
as  a function  of  Ng  = Ng  - Ng.  The  quantity  Ng  is  usually 
called  the  signal  excess,  although  it  is  defined  in  terms  of 
signal  and  noise  level  differences.  If  the  dependence  of  Ng 
on  record  duration  is  known,  then  one  curve  will  likely  suffice 
for  all  of  the  duration  options,  provided  that  no  other  display 
parameters  are  changed. 

A single  transition  curve  does  not  completely  characterize 
the  combined  operation  of  the  components  to  the  right  of  the 
dashed  line.  A different  curve  could  obtain  if  the  power  fluc- 
tuations of  the  inputs  were  distributed  differently,  or  if  the 
time  scales  of  those  fluctuations  were  different,  as  will  be 
seen  in  Section  7.0 
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It  is  difficult  to  obtain  data  for  a transition  curve 
under  actual  operating  conditions  because  of  the  problem  of 
determining  the  signal  level  at  the  output  of  the  receiver 
bandpass  filter  when  it  is  much  less  than  the  noise  level. 
Figure  1.2  shows  a means  for  obtaining  the  curve  under 
laboratory  conditions.  'Signal  levels  would  be  varied,  and 
numerous  observations  would  be  made  by  several  or  more  opera- 
tors. The  relevance  of  the  curve  to  a particular  environment 
will  depend  on  the  degree  to  which  the  statistics  of  the  inputs 

approach  those  under  operating  conditions. 

* 

Given  the  appropriate  transition  curve  as  a function  of 
Ng,  the  remaining  task  is  to  find  the  value  of  NE  achieved  by 
the  sonar  system  under  specified  operating  conditions.  This 
is  usually  done  by  means  of  the  sonar  equation 


Nr 


LS  " Nw  - LN  + na  - ND 


(1.1) 


where  Lg  is  the  source  level  in  the  receiver  passband 
Nw  is  the  propagation  loss 

LN  is  the  noise  level  at  the  array  in  the  receiver 
passband 

Na  * 10  log  Gg  - 10  log  Gn 


Gg  is  the  average  signal  power  gain  of  the  array  over 
the  receiver  passband 


'N 


is  the  average  noise  power  gain  of  the  array  over 
the  receiver  passband 
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If  Ng  is  identically  zero  in  (1.1),  the  result  is  the  sonar 
detection  threshold  equation;  solving  that  equation  for  Nw 
gives  the  figure  of  merit,  the  value  of  Nw  for  which  the 
probability  of  detection  is  that  specified  for  Np.  Using 
this  particular  value  of  Nw  and  the  appropriate  transmission 
loss  curve  gives  the  detection  threshold  range. 

If  the  sonar  system  is  operated  in  different  environments, 
then  some  of  the  sonar  equation  values  on  the  right  side  of  (1.1) 
are  not  known  a priori,  and  they  can  be  regarded  as  random 
variables.  In  that  case,  it  is  clear  from  (1.1)  that  N£  is 
also  a random  variable.  In  that  case,  the  single-look  prob- 
ability of  detection  is  given  by 


where 


Pn  = E{g(N„)}  = / dx  g(x)f„(x) 


E is  the  expectation  operator 


g(  ) is  the  function  specifying  the  transision 
curve 


fE(  ) is  the  probability  density  function  for  Ng 

If  Ng  is  broadly  distributed  compared  to  the  spread  of  the 
transision  curve,  then 


(1.2) 


pD  - / dx  fg(x) 


(1.3) 
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This  is  equivalent  to  assuming  that  g(x)  is  a unit  step  function 
u(x).  If  the  terms  on  the  right  of  (1.1)  are  statistically  inde- 
pendent, then  the  distribution  of  N£  is  readily  determined  if 
the  distributions  of  the  terms  are  known.  However,  the  terms 
may  not  be  independent;  for  example,  Nw  and  may  be  dependent 
to  the  extent  that  the  noise  sources  are  in  the  propagation 
medium. 


Transition  curves  can  also  be  obtained  by  employing  the 
methods  of  statistical  detection  theory.  This  requires  that 
a circuit  analog  be  developed  for  the  components  to  the  right 
of  the  dashed  line  in  Figure  1.1.  The  classical  single-channel 
analog  (II. B,  Reference  1)  includes  a square-law  rectifier, 
an  averager,  and  a threshold  detector.  A multi-channel  analog, 
described  in  Section  2.0,  is  applicable  to  the  analysis  of  a 
wider  variety  of  input  conditions. 

The  application  of  detection  theory  also  requires  charac- 
terizing the  inputs  statistically.  In  Section  6.1,  the  transition 
curve  is  obtained  for  a case  in  which  the  Inputs  are  stationary 
normal  stochastic  processes.  In  Section  6.2,  the  inputs  have 
the  form  V G(t),  where  P is  a non-negative  random  variable  with 
the  dimensions  of  power,  and  G(t)  is  a stationary  normal  random 
process  with  zero  mean  and  unit  variance.  This  analysis  leads 
to  a calculation  equivalent  to  (1.2).  Several  examples  are 
investigated  in  which  P is  a random  variable  with  a gamma  dis- 
tribution. For  these  cases,  the  integral  in  (1.3)  yielded 
convenient  closed  forms;  the  accuracy  of  this  approximation 
was  investigated  and  found  to  be  good. 
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A more  general  class  of  Inputs  has  the  form  ^ P(t )G(t ) , 
where  P(t),  called  a power  envelope  process.  Is  a non-negative 
stationary  random  process  with  dimensions  of  power.  In  Section 
4.2,  it  is  shown  that  the  statistical  properties  of  this  class 
agree  with  the  results  of  careful  measurements  of  ambient  noise. 
Section  4.1  discusses  the  relevance  of  stationary  processes  to 
cases  in  which  the  inputs  have  cyclical  non-stationarities 
whose  periods  are  very  large  compared  to  the  post-rectification 
averaging  time.  Section  4.3  defines  and  analyzes  two  classes 
of  power  envelope  processes. 

Section  7-0  produces  transition  curves  for  cases  in  which 
the  noise  inputs  have  power  envelopes  that  are  fully  correlated. 
For  performance  predictions,  these  curves  would  be  utilized  by 
employing  (1.1)  to  calculate  the  value  of  Ng.  In  this  case, 

ND  is  that  pertaining  to  stationary  normal  inputs;  it  is  cal- 
culated from 

Nd  * 10  log  y - 5 log  WT  (1.4) 


where  N(-y)  - Pp 

N(  ) is  the  normal  probability  distribution 
function  for  a zero-mean,  unit-variance, 
random  variable 

Pp  is  the  probability  of  false  alarm  with 
normal  noise  inputs 

W is  the  noise  equivalent  bandwidth,  defined 
below  (5.17) 

T is  post-rectification  averaging  time 
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At  the  present  time,  (1.2)  represents  the  state-of-the-art 
of  dealing  with  the  effects  of  variability  of  sonar  parameters, 
alghough  (1.3)  is  more  frequently  employed.  Either  applies  if 
one  or  more  of  the  parameters  is  a random  variable,  or  if  the 
relaxation  time  is  long  compared  to  the  post-rectification 
averaging  time.  The  method  developed  in  this  report  does  not 
assume  those  conditions,  and  is  capable  of  showing  the  depen- 
dence of  performance  on  the  relaxation  time  of  the  power 
envelopes.  Figure  7-3,  for  example,  shows  transition  curves 
for  different  values  of  T * D,  where  D is  the  relaxation  time 
of  the  noise  power  envelope. 

The  remainder  of  the  report  develops  the  relations  required 
for  calculating  the  transition  curves  for  the  case  of  the  fluc- 
tuating power  envelope.  Section  3*0  develops  a means  for 
approximating  the  probability  density  of  the  averager  output 
based  on  a polynomial  expansion  and  a basis  density  function. 
Section  5-0  derives  low-order  moments  of  the  averager  output; 
these  moments  are  required  for  the  calculation  of  the  coefficients 
of  the  polynomial  expansion.  An  alternative  to  this  approach 
is  discussed  in  the  Appendix. 
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2.0  MULTICHANNEL  ANALOG 

2.1  Introduction 

Certain  functions  performed  by  sonar  systems  require 
simultaneous  processing  of  multiple  inputs.  One  example  is 
the  processing  of  multiple  broadband  preformed  beam  signals 
to  be  displayed  in  a beam-time  format.  Another  example  is  a 
frequency  analyzer  utilizing  either  a bank  of  fixed  filters 
or  more  likely  a discrete  Fourier  transform  algorithm.  Here 
the  output  format  is  bin  or  filter  number  (representing  a 
frequency  band)  versus  time.  In  these  examples,  the  channels 
are  discrete;  however,  each  has  a continuous  analog  in  which 
either  a beam  is  swept  in  azimuth,  or  an  oscillator  is  swept 
in  frequency  to  heterodyne  signals  into  a single  fixed  pass- 
band.  These  continuous  versions  can  be  represented  in  discrete 
forms  with  contiguous  multiple  channels  with  the  appropriate 
bearing  or  frequency  resolution. 

The  multichannel  analog  was  initially  employed  [Ref.  1, 
Section  II]  to  predict  the  performance  of  an  operator  using  a 
multichannel  display.  It  was  applied  to  the  case  of  a broad- 
band multibeam  sonar  with  a beam-time  display  format.  With 
this  approach,  it  was  possible  to  predict  the  effects  of  acoustic 
interference  produced  by  nearby  ships  on  sonar  detection  perform- 
ance. The  limiting  case  of  infinite  averaging  time  has  been 
labeled  the  bias-limited  case  [2] . A more  compact  derivation 
of  the  performance  of  the  multichannel  analog  is  given  in 
Ref.  3. 


A more  detailed  analysis  of  the  detection  performance 
of  an  operator  using  Lofar  [4]  was  based  on  certain  results 
of  neurophysiological  and  psychophysical  studies  of  the  visual 
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perception  system.  One  result  of  that  analysis  showed  that  the 
effective  frequency  resolution  of  the  operator's  receptive* 
fields  could  be  larger  than  that  of  the  Lofar  processor,  and 
that  detection  performance  could  be  poorer  than  that  expected 
on  the  basis  of  the  analyzer  bandwidth.  If  the  multichannel 
analog  is  employed,  a first-order  correction  of  this  discre- 
pancy can  be  made  by  selecting  the  larger  value  of  bandwidth. 

A schematic  of  the  multichannel  analog  is  shown  in 
Figure  2.1.  The  representation  is  vectorial  in  that  it 
represents  multiple  signal  flow  and  multiple  (but  similar) 
components  in  single  blocks.  The  outputs  of  the  averagers 
are  a set  of  quantities  Y = {Y^.  In  the  classical  circuit 
analog  [Ref.  1,  pp.  8-14],  each  of  these  would  be  compared  to 
a fixed  threshold  to  form  the  basis  for  a detection  decision. 

In  the  multichannel  analog,  however,  a quantity  Z^  is  formed 
for  each  channel  that  is  the  difference  between  the  averager 
output  Y^  and  a quantity  that  is  a linear  combination  of 
the  averager  outputs  Y^  i ? k.  If  the  quantity  Z^  Is  greater 
than  zero,  a detection  would  be  declared.  This  scheme  is  equi- 
valent to  that  described  in  Refs.  1 and  3 in  which  Y^  is  compared 
to  a threshold  Wk.  The  advantage  of  the  formulation  shown  in 
Figure  2.1  is  that  the  basis  for  decisions  is  a single  random 
vector  Z instead  of  the  pair  Y and  W. 

2.2  Detection  Analysis 

If  the  post-rectification  averaging  is  uniform  for  a 
time  period  T,  then  the  output  of  the  1th  averager  can  be  ex- 
pressed as 

*The  receptive  field  of  a neuron  in  the  visual  system  is  defined 
as  that  portion  of  the  retina  (or  the  corresponding  visual 
field)  that  produces  a response  in  the  neuron  when  stimulated 
by  light. 
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Y,  (t)  = T"1  / du  X,2(u) 

1 t-T  1 


(2.1) 


If  the  Inputs  are  stationary,  then  the  steady-state  probability 
measures  are  Invariant  with  time,  and  it  will  be  convenient  to 
evaluate  the  output  at  time  T: 


Y,  = T"1  I du  X2(u) 
1 0 1 


(2.2) 


If  no  signal  is  present  in  channel  i,  then  the  output  to 
its  squarer  is 


Xj_(t)  * Ni(t) 


and  the  output  of  its  averager  is 


(2.3) 


Y, (t)  = T”1  I du  N2 (U) 
1 0 1 


(2.4) 


And  if  a signal  is  present  in  channel  0 (zero),  then  the  input 
to  its  squarer  is 


XQ(t)  - NQ(t)  + S(t) 


(2.5) 


and  the  output  of  its  averager  is 


YQ(t)  - T-1  / du  [Nq (u)  + 2NQ(u)S(u)  + S2(u)] 


(2.6) 


where  S(t)  is  the  signal  component  of  the  input  to  channel  0. 
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Per  the  discussion  in  Section  2.1,  the  test  statistic 
for  channel  0 is 


1=1  1 1 


(2.7) 


(2.8) 


where  c^  is  the  weighting  coefficient  for  channel  i.  To 
simplify  the  notation,  the  subscript  zero  will  be  suppressed 
except  to  designate  the  noise  component  of  the  input  to  that 
channel.  Substituting  (2.4)  and  (2.6)  in  (2.8)  gives 


T n 

Z = T-1  / du  [S2 (u)  + 2S(u)N_(U)  - Z c,N*(u)] 
0 u 1=0  1 1 


(2.9) 


where  the  noise  term  for  the  channel  of  interest  is  included  in 
the  summation  by  setting  cQ  = -1. 

The  probability  that  the  test  statistic  Z exceeds  the 
(zero)  threshold  is 


P(Z>0)  = / dzf7(z) 


(2.10) 


where  fz(  ) is  the  probability  density  function  for  the  test 
statistic  Z.  It  will  be  convenient  to  consider  the  standardized 
random  variable 


(2.11) 
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where  m £ and  are  the  mean  value  and  standard  deviation  of  Z 
respectively.  Then  the  probability  that  the  test  statistic  Z 
exceeds  the  threshold  is 


P(Z>0)=  P(Q  > - mz  t oz) 


(2.12) 


where  fq(q)  = + mZ^  is  probability  density  function 

for  the  standardized  random  variable  Q.  The  (conditional)  prob- 
ability of  detection  is  then 

PD  = P(Z>0|ps>0 1 


mJ  dqfQ(q|p^>0)  (2.13) 

“mzs*azs 


where  pg  is  the  mean  square  value  of  S(t) 
mzs  * E(Z|ps>0) 

°zs  * »(zIps>0) 
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I 

And  the  probability  of  false  alarm  is 

Pp  = P(Z>0|PS  = 0) 

00 

= J dqfQ(q|Ps  = 0)  (2.14) 

_mZN*aZN 

where  MZN  = E(Z|ps  = 0) 

°2N  ' 0(ZIPS  * °> 

The  bases  for  further  analysis  are  (2.13)  for  the  probability 
of  detection  on  a single  observation,  (2.14)  for  the  probability 
of  false  alarm  on  a single  observation,  (2.11)  for  the  standardized 
random  variable,  and  (2.9)  for  the  test  statistic. 
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3.0  OUTPUT  DENSITY  FUNCTION 


As  shown  in  the  previous  section,  the  calculation  of 
the  probabilities  of  detection  and  false  alarm  requires  the 
probability  distribution  of  the  standardized  form  of  the  test 
statistic.  The  latter  is  a linear  combination  of  the  outputs 
of  Integrators  with  non-Gaussian  inputs.  The  general  problem 
of  determining  the  required  distribution  function  is  difficult. 
Papoulis  (Ref.  5>  P*  324)  states: 


"The  determination  of  the  distribution  Fo(s) 
of  S [integrator  output]  is,  in  general,  hopelessly 
complicated.  After  all,  S is  the  limit  of  a sum  of 
[random  variables],  and  as  we  know,  even  to  find  the 
distribution  of  the  sum  of  only  two  random  variables 
is  not  a trivial  matter.  For  this  reason  we  shall 
not  be  concerned  with  FgCs),  but  shall  determine  only 
the  mean  and  variance  of  S." 


The  same  remark  applies  as  well  to  the  probability  density 
function  of  S. 


Cramer  has  shown  (Ref.  6,  Ch.  12,  Sec.  6)  that  a prob- 
ability density  function  can  be  determined  from  its  moments. 

A density  function  g(x)  can  be  expanded  in  the  form 

00 

g ( x ) » f (x)  E b .p. (x)  (3-1) 

i*0  1 1 

where  f(x)  is  a selected  density  function 

bj^  are  constants  to  be  determined 

P1(x)  is  an  lth-order  polynomial  satisfying  the 
condition 
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00 

f dx  f (x)pm(x)pn(x)  = 
—00 


!1  for  m = n 
0 for  m / n 


(3-2) 


The  coefficients  of  the  polynomial  are  determined  by  the  moments 
of  f(x),  and  the  coefficients  explicit  in  (3-1)  are  given  by 


n 


-/ 


00 

dX 


Pn  (x)g(x) 


(3.3) 


Since  Pn(x)  is  a polynomial  in  x,  it  is  seen  that  bn  is  determined 
Jointly  by  the  coefficients  of  the  polynomial  and  the  moments  of 
g(x). 


An  approximation  of  g(x)  is  obtained  by  using  a finite 
number  of  terms  in  the  expansion.  One  expects  that  the  approxi- 
mation will  be  good  if  f(x)  is  itself  a reasonably  good  approxi- 
mation to  g(x).  If  the  series  (3.1)  is  terminated  at  a finite 
number  of  terms,  then  the  coefficients  of  like  powers  of  x can 
be  summed  to  obtain  an  expansion  of  the  form 

h(x)  = f(x)P(x)  (3.4) 

n i 

where  P(x)  = E a.  x 
i*0  1 

Since  h(x)  is  to  approximate  a probability  density  function,  it 
appears  reasonable  to  require  that 


f h(x)xl 


'dx  = a^,  k * 0,  1,  2,  ...  n 


(3.5) 


r 
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where  ak  is  the  kth  moment  of  g(x).  Equations  (3.5)  require 
that  / ■ dxh(x)  = 1,  one  requisite  of  any  probability  density 
function,  and  that  the  first  n moments  of  h(x)  equal  those  of 
g(x).  Substituting  (3-*0  in  (3-5)  and  evaluating  the  integrals 
gives 

n 

£ ^i  ^k+  i — ak  (3*6) 

where  m^  is  the  Jth  moment  of  the  basis  function  f(x),  and  aQ  = 1. 

In  Section  2.0,  the  probability  of  detection  is  expressed  as 
an  Integral  of  the  probability  density  function  of  a standardized 
random  variable;  for  such  a variable,  the  mean  is  zero,  the  vari- 
ance one,  and  all  moments  are  central.  If  both  g(x)  and  the  basis 
function  f(x)  are  densities  for  standardized  random  variables,  then 
the  matrix  and  vector  representing  (3-6)  for  n = 3 are  as  shown 
below: 


1 

0 

rH 

m3 

1 

0 

1 

m3 

m^ 

0 

1 

m ^ 

m5 

1 

m3 

m4 

m5 

m6 

a3 

It  is  easy  to  show  that  the  coefficients  a^,  a2>  and  a^  are 
proportional  to  the  difference  of  the  third  moments  - m y For 
example,  to  solve  for  a2  by  Cramer's  rule,  replace  the  third  column 
of  the  matrix  by  the  vector,  and  then  subtract  the  elements  of  the 
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first  column  from  the  new  third  column.  The  result  is  a column 

with  all  zeros  except  for  the  fourth  element  - m^.  If  the 

determinant  is  expanded  in  terms  of  the  elements  of  this  column  i 

and  their  minors,  it  is  apparent  that  a2  is  proportional  to  the 

difference  of  third  moments  - m^.  A similar  procedure  applies 

to  a^  and  a^.  ’ 

The  method  is  attractive  in  that  it  requires  relatively  few 
moments  of  the  test  statistic.  Its  mean  and  variance  are  required 
for  (2.13)  and  (2.14)  for  the  probabilities  of  detection  and 
false  alarm.  For  the  case  just  discussed,  the  third  moment  is 
also  required.  Although  the  mean  value  is  easily  obtained,  the 
derivation  of  the  variance  and  higher  moments  requires  considerable 
effort.  For  the  basis  function  f(x),  on  the  other  hand,  the 
number  of  moments  required  is  twice  that  for  g(x).  However,  the 
. problem  of  obtaining  these  moments  is  usually  not  difficult. 

For  many  density  functions,  the  characteristic  function  is  avail- 
able for  generating  the  moments,  or  in  some  cases,  formulas  for 
the  moments  are  given,  as  in  P.cference  7.  In  any  case,  (3*5) 
can  be  applied  with  h(x)  replaced  by  f(x),  and  by  n^. 

For  n * 3,  the  polynomial  may  be  expressed  as 

_i  3 t 

P(x)  - 1 + (o  - m ) A E fi  x1  (3.8) 

J J i=0  1 

where  A * fi  + figin^ 
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- m*  - 1 


m5  - m3mi} 


- m3 


mjj  “ 1 


+ 1) 


m6  - 


ra5  — 


m.-  - m_ 
o 3 


m._  - m., 
o 3 


mi|  ” 1 


In  many  cases  of  interest  g(x)  will  be  approximately  normal, 
and  the  normal  density  function  is  an  appropriate  basis  for  the 
expansion.  For  this  case,  all  odd-order  moments  of  f(x)  are  zero, 
m^  * 3,  and  mg  ■ 15.  Evaluating  (3.8)  and  substituting  the  result 
in  (3.1*)  gives 


h(x)  * n(x)  II  - a^(x/2  - x3/6)] 


(3.9) 


where  n(x)  is  the  density  for  a normal  variate  with  zero  mean  and 
unit  variance.  The  coefficients  a.^  through  a^  are  seen  to  be  pro- 
portional to  (a^  - m^),  since  m^  * 0 for  this  basis  function. 
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In  view  of  the  relationship  between  the  Hermite  poly- 
nomials and  the  derivatives  of  n(x),  it  can  be  shown  that 
an  alternative  form  for  (3*9)  is  [Ref.  6,  pp.  222-223] 


i 

» 


h(x) 


(3.10) 


Substituting  this  result  in  (2.12)  and  evaluating  the  integrals 
gives 


P(Z>0) 


mZ 

N(— =■) 
°Z 


y3Z 

3 

°z 


(3-11) 


i 

where  N(  ) is  the  probability  distribution  function  for  a normal 
random  variate  with  zero  mean  and  unit  variance. 


This  is  a very  convenient  form  since  tables  of  both  N(  ) and 
n(i)(  ) are  available. 

In  Section  2.0,  the  output  of  the  multichannel  analog  (2.7) 
was  expressed  as 


Z = Y - W (3.12) 

where  Y is  the  output  of  the  averager  of  the  channel  of  interest 

W is  a threshold  function  synthesized  from  the  averager 
outputs  of  other  channels. 

The  former  (Y)  is  non-negative,  and  the  latter  (W)  will  also  be 
if  none  of  the  weighting  coefficients  aj^  in  (2.8)  are  negative. 
Given  these  conditions,  the  probability  density  of  Z can  be 
expressed  as 


i 

■ 

\ 
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00 

fz(z)  = f dxfyw(x  + z,  x),  z >0  (3.13) 

t 

00  , 

= / dxfyw(x,  x - z),  z<0  (3-1*0 

where  fyw  ( , ) Is  the  joint  density  function  for  Y and  W.  Unless 
Y and  W are  completely  dependent,  the  variate  Z will  be  bipolar; 
furthermore,  the  density  function  is  continuous  at  z = 0 as  can 
be  seen  from  (3*13)  and  (3.1*1). 


To  illustrate  how  the  density  function  of  the  test  statistic 
Z compares  to  that  of  the  averager  output  Y of  the  channel  of 
interest,  the  following  extreme  case  is  considered.  The  averager 
output  Y is  chi-square  with  two  degrees  of  freedom,  and  the 
threshold  function  W is  one-half  of  a chi-square  variate  with 
four  degrees  of  freedom,  the  result  of  averaging  the  outputs  of 
the  averagers  of  two  adjacent  channels.  These  processes, 
assumed  independent,  would  result  with  an  extremely  small  time- 
bandwidth  product.  Applying  (3-13)  and  (3.1**)  to  this  case 
yields  the  density  function  plotted  in  the  lower  part  of 
Figure  3.1.  The  function  is  continuous  everywhere,  as  is  its 
first  derivative.  The  upper  part  of  Fig.  3.1  shows  the  density 
functions  for  Y and  -W.  The  former  is  discontinuous  at  zero, 
as  is  the  first  derivative  of  the  latter.  The  result  of  the 
combination  via  (3*12)  yields  a test  statistic  that  is  more 
nearly  normal  than  the  output  of  the  averager  of  the  channel 
of  interest.  For  the  example  selected  coefficients  of  skewness 
and  kurtosis  are  0.82  and  6 respectively  for  the  former,  and 
1 and  9 respectively  for  the  latter. 
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x,  distribution  parameter  of  Y and  -W. 
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z,  distribution  parameter  of  Z. 


FIGURE  3.1  Probability  Density  Functions. 
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4.0  A COMPOUND  RANDOM  PROCESS 

4.1  Introduction 

The  results  of  several  experimental  studies  (Refs.  8 
and  9)  of  ambient  noise  may  be  summarized  as  follows: 

1.  Estimates  of  first-order  probability  measures 
obtained  from  relatively  short  segments  (com- 
pared to  certain  relaxation  times  of  the  process) 
of  sample  functions  supported,  more  often  than 
not,  the  hypothesis  that  the  amplitude  of  noise 
in  a specified  bandpass  is  distributed  normally. 

2.  The  number  of  departures  from  first-order 
normality  increased  with  sample  length  (Ref.  8, 

Tables  V and  VI.) 

3.  Estimates  exhibited  more  variability  than 
consistent  with  a stationary  Gaussian  pro- 
cess. 

The  first-order  probability  measures  include  the  distribution 
function  and  several  of  its  low-order  moments. 

It  will  be  shown  in  the  next  section  that  similar  results 
could  be  expected  by  observing  a stationary  compound  random 
process  that  is  conditionally  non-stationary  normal. 

A complete  process  model  for  ambient  noise  would  be  non- 
stationary  because  of  dependence  on  phenomena  with  cycles  such 
as  diurnal,  lunar,  or  annual.  Figure  4.1  shows,  for  example, 


26 


RELATIVE  LEVEL  WIND  SPEED 


titikU&Utsir&ift.  - 


Bolt  Beranek  and  Newman  Inc. 


Report  No.  3656 


I 

I 

I 

I 

I 


the  dependence  of  noise  power  level  on  time  of  day.  There  may 
also  be  long-term  trends  due  to  evolution  of  the  generating 
phenomena,  for  example,  the  active  ship  population. 

If,  however,  the  variation  of  process  parameters  is  small 
during  a period  of  interest,  the  process  may  be  considered  as 
locally  stationary  during  that  period,  and  approximated  by  a 
stationary  process  with  parameters  that  are  appropriate  to  the 
time.  If  the  shortest  cycle  is  diurnal,  and  if  the  post- 
rectification averaging  period  does  not  exceed  an  hour,  this 
approach  should  be  entirely  satisfactory  for  the  objectives  of 
this  investigation. 

A further  requirement  for  the  model  is  that  its  mean  value 
be  zero.  This  requirement  would  be  met  by  any  stationary  process 
that  is  the  output  of  either  a high-pass  or  band-pass  network 
or  operation.  For  example,  the  rth  complex  coefficient  of  the 
discrete  Fourier  transform  is 

n-1 

A * Z X.  exp(-  J2irrkn“l)  (4.1) 

r k-0  K 


where  Xk  is  the  kth  sample  of  the  input  (assumed  stationary)  in 
the  frame.  The  expected  value  of  the  coefficient  is 

n-l 

E( A ) « mY  E exp  (-  J2irrkn_1)  (4.2) 

r * k-0 

where  m^  - E(Xk),  all  k.  It  is  easily  seen  from  (4.2)  that  for 
r - mn,  n - 0,  1,  2,  3,  ....,  that  the  value  of  the  indicated 
sum  is  n.  This  set  of  outputs  represents  the  low-pass  output 
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and  its  aliases.  For  other  values  of  the  index  r,  the  geometric 
progression  defined  by  the  sum  is 


l = 1 ~ exp  (-  J 2irr) 

1 - exp  (-  j2Trrn_1) 


(4.3) 


The  numerator  is  zero  for  all  r,  and  the  denominator  is  not  ex- 
cept for  r * mn.  Thus 

E ( Ap ) = nmx , r=mn,  n = 0,  1,  2,  3,  .... 

= 0,  r an  integer  / mn  (4.4) 

This  set  of  outputs  represents  a set  of  bandpass  outputs  and 
their  aliases.  This  result  has  been  demonstrated  experimentally 
(Ref.  9,  Page  68). 

4.2  A Class  of  Compound  Processes 

A fairly  broad  class  of  random  processes  exhibiting  the 
behavior  described  in  the  previous  section  is  the  compound  pro- 
cess 


N(t ) =^P(t)  G(t ) 


(4.5) 


where  P(t)  is  a non-negative  stationary  random  process. 

G(t ) is  a zero-mean,  unit-variance,  Gaussian  pro- 
cess, statistically  independent  of  P(t). 

The  class,  briefly  described  in  Pages  153-154  of  Ref.  10,  is 
broad  in  that  the  only  constraint  on  P(t)  is  that  it  not  be 
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negative.  It  is  clear  from  (4.5)  that  given  P(t),  the  process 
N(t)  is  conditionally  non-stationary  Gaussian. 

The  compound  process  defined  by  (4.5)  includes  two  sub- 
cases of  more  than  passing  interest.  In  one  of  these,  the 
power  envelope  P(t)  is  a random  variable  P independent  of  time; 
in  this  case,  the  compound  process  is  non-ergodic,  since  time 
averages  of  individual  members  would  generally  not  be  equal  to 
ensemble  averages.  Since ^"p~ is  not  negative,  the  compound  pro- 
cess ■^T’G(t)  is  a special  case  of  a spherically  invariant  process 
A G(t ) , in  which  A is  a random  variable  that  Is  not  necessarily 
non-negative  (Ref.  11).  In  the  second  case,  the  power  envelope 
is  a deterministic  constant,  and  the  compound  process  degenerates 
to  a Gaussian  process  (ergodic,  and  non-compound).  Both  of  these 
special  cases  will  be  considered  further  in  Section  6.0. 

The  probability  density  function  for  N(t)  can  be  expressed 
as 


where  fp(  ) is  the  probability  density  function  for  P(t).  It  is 
clear  from  (4.6)  that  the  function  has  even  symmetry  around  n = 0; 
hence,  all  odd-order  moments  about  zero  are  zero  including  the 
first,  which  satisfies  one  of  the  requirements  stated  in  the  pre- 
vious section.  Furthermore,  all  even-order  moments  are  moments 
about  the  mean. 

For  purposes  of  illustration,  consider  a case  in  which  P(t) 
is  a chi-square  variable  with  four  degrees  of  freedom.  Applica- 
tion of  (4.6)  to  this  case  yields 


30 


Bolt  Beranek  and  Newman  Inc. 


Report  No.  3656 


fN(n)  = \ (1  + | n | ) exp  ( - | n | ) 


(4.7) 


which  Is  plotted  In  Figure  4.2  along  with  the  density  for  a 
normal  variate  with  the  same  variance.  The  departure  from 
normality  Is  quite  evident  from  the  density  functions;  the 
departure  would  be  less  evident  from  the  corresponding  distri- 
bution functions. 

Even-order  moments  are  readily  calculated  from  the  corres- 
ponding powers  of  (4.5),  the  square  of  which  is 


N2 (t ) * P(t )G2 (t ) 


(4.8) 


If  the  relaxation  time  of  P(t)  is  very  large  compared  with  that 
of  the  normal  process  G(t),  then  it  seems  reasonable  to  regard 
the  former  as  a power  envelope  process.  The  second  moment  of 
the  compound  process  is 


E [N2 (t ) ] = E [P(t ) ] E [G2 (t ) ] 


(4.9) 


where  p = E[P(t)] 


The  fourth  central  moment  of  N(t)  is 


- E [G4 (t ) ] E [P2 (t ) ] 


3(p2  + a2) 


where  a2  is  the  variance  of  P(t). 


(4.10) 


compound  process 
normal  process 
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The  departure  from  first-order  normality  of  the  compound 
process  can  be  gauged  by  comparing  its  moments  to  those  of 
corresponding  order  of  a normal  variate.  For  a zero-mean  normal 
variate,  all  odd-order  moments  are  zero,  as  is  the  case  for  the 
compound  process  given  by  (4.5);  thus,  for  both  of  these  processes 
the  coefficient  of  skewness  * a3  is  zero.  Figure  4.3(a)  shows 
the  distribution  of  sample  estimates  of  that  coefficient  obtained 
by  Arase  and  Arase  for  deep-sea  ambient  noise.  The  distribution 
is  fairly  even  about  zero,  and  about  12  percent  of  the  sample 
values  fill  outside  of  the  95  percent  confidence  limits.  For  a 
normal  variate,  only  5 percent  of  the  values  are  expected  outside 
these  limits. 

The  kurtosis  of  a distribution  is  defined  as  p4  * a1*. 

Figure  4.3(b)  shows  the  distribution  of  kurtosis  estimates  for 
deep-sea  ambient  noise  obtained  by  Arase  and  Arase.  In  this 
case  about  28  percent  of  the  values  are  outside  the  95  percent 
confidence  limits,  whereas  only  5 percent  would  be  expected  for 
a normal  variate.  Utilization  of  (4.9)  and  (4.10)  shows  that 
the  kurtosis  for  the  compound  random  process  is 

P4  * orj  - 3 [1  + (a/p)2]  ' (4.11) 

Thus,  if  the  Variance  of  the  power  envelope  process  is  greater 
than  zero,  the  kurtosis  of  the  compound  process  is  greater  than 
three,  the  kurtosis  of  a normal  population.  For  the  compound 
process  whose  probability  density  function  is  given  by  (4.7), 
the  value  of  kurtosis  is  (4.5)  The  kurtosis  estimates  discussed 
above  were  obtained  from  data  periods  of  10,  20,  and  40  seconds, 
and  the  normal  hypothesis  was  rejected  more  often  with  the 
longer  data  periods.  It  seems  reasonable  to  expect  that  larger 
estimates  would  be  obtained  with  longer  sample  periods.  These 


33 


Bolt  Beranek  and  Newman  Inc. 


Report  No.  3656 


» 

' f 


! 

i: 

i 

i 

i 

i 

i 

i 

i 

i 


experiments  were  conducted  with  single  hydrophones;  beam  noise 
data  might  be  expected  to  exhibit  greater  variability  because 
of  the  highly  uneven  weighting  of  responses  to  the  noise  sources, 
which  has  the  effect  of  reducing  the  size  of  the  noise  source 
population. 

Certain  of  the  temporal  characteristics  of  the  compound 
random  process  are  revealed  by  low-order  moment  functions, 
easily  derivable  from  (4.5).  The  autocovariance  function  for 
the  process  is 


Kn(u-v)  = E[N(u)N(v)] 

= E[^P(u)P(v)]  p(u-v)  (4.12) 

where  p(  ) is  the  autocovariance  coefficient  function  for  the 
normal  process  G(t).  If  the  relaxation  time  of  the  P-process  is 
much  longer  than  that  of  the  G-process,  then 

Kn(u-v)  » E [P(u) ] p(u-v)  = pp(u-v);  (4.13) 

that  is,  the  G-process  essentially  determines  the  relaxation  time 
of  the  compound  process  and  thereby  its  bandwidth.  It  is  usually 
this  bandwidth  that  serves  as  the  basis  for  decisions  regarding 
sample  intervals  and  data  periods  for  experimental  efforts. 

The  data  periods  based  on  this  bandwidth  may  not  be  very  long 
compared  to  the  longer  relaxation  times,  in  which  case  the  number 
of  statistically  independent  samples  will  not  be  adequate  to 
obtain  good  parameter  estimates. 
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Since  the  receiver  analog  includes  squarers,  the  properties 
of  Q(t)  = N2(t)  are  important;  the  second  moment  function  for 
the  square  of  N(t)  is 

Mq(u-v)  = E [N2 (u)N2 (v) ] 

= E [P(u)P( v) ] E [G2 (u)G2 ( v) ] (4.14) 

This  moment  function  is  seen  to  involve  a fourth-order  moment 
function  for  N(t);  however,  the  form  is  degenerate  in  that  there 
is  only  one  time  difference,  u-v.  The  second  factor  of  (4.14) 
is  a fourth-order  moment  function  (degenerate)  for  a normal  pro- 
cess. This  function  can  be  expressed  in  terms  of  lower-order 
properties  by  using  (7-28)  of  Ref.  12: 

Mq(u-v)  = (p2  + K(u-v) ] [1  + 2p2(u-v)]  (4.15) 

where  K(u-v)  is  the  autocovariance  function  for  the  power  en- 
velope process.  Given  (4.15)  and  (4.9),.  it  follows  that  the 
autocovariance  function  for  the  square  of  the  compound  process 
is 


Kq(u-v)  = K(u-v)  + 2 [p2  + K(u-v) ] p2 (u-v)  (4.16) 

Given  the  stipulation  on  relaxation  times, 

Kq(u-v)  = K(u-v)  + 2 (p 2 + a2)p2(u-v)  (4.17) 

Given  that  the  variance  a2  of  the  power  envelope  process  is  not 
negligible  compared  to  p2,  then  the  relaxation  time  of  the  square 
of  the  compound  process  is  essentially  determined  by  that  of  the 
power  envelope  process.  This  is  just  the  reverse  of  the  situa- 
tion pertaining  to  the  first  power  of  the  process,  and  the 
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ramifications  with  regard  to  the  effectiveness  of  post-rectification 
smoothing  for  either  detection  or  power  estimation  are  considerable. 

Figure  4.4  shows  a time  record  of  the  square  root  of  suc- 
cessive estimates  of  the  mean  squared  amplitude.  There  is  an 
overall  trend  in  data  correlated  with  a reduction  of  local  wind 
speed.  Within  this  period  about  three  cycles  of  shorter  periods 
are  evident.  (These  periods  are  much  longer  than  the  data  averag- 
ing period,  in  this  case  20  seconds.)  Similar  behavior  would  be 
exhibited  by  samples  from  a compound  process  with  a power  envelope 
process  exhibiting  two  relaxation  times,  the  longer  of  which  is 
related  to  the  relaxation  time  of  wind  speed. 

All  of  the  data  cited  lend  support  to  the  consideration  of 
(4.5)  as  a reasonable  model  for  ocean  ambient  noise. 

t 

Analysis  of  the  performance  of  the  multichannel  analog  requires 
specification  of  a vector  random  process,  that  is,  of  a set  of 
random  processes.  The  set  is  described  by 


Ni(t)  = G1(t),  i - 1,  2,  3,  .....  n (4.18) 

where  P1(t)  is  a non-negative  stationary  random  process  with 
average  value  p^ 

G^t)  is  a zero-mean,  unit-variance,  Gaussian  process, 
statistically  independent  of  P1(t),  all  i,  and 
statistically  independent  of  Gj(t),  all  J except  i. 
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The  assumed  mutual  Independence  of  the  Gaussian  processes  would 
be  approximated  by  beamformer  outputs  In  a band  of  frequencies 
near  the  design  frequency  of  the  array.  Note  that  the  power 
envelope  processes  P^Ct)  may  be  mutually  dependent. 

The  mean  value  of  the  1th  member  is 

E[Nj_(t)]  = 0 (4.19) 


and  the  cross-covariance  is 

Kij(u-V)  “ P^^u-v),  J = 1 

- 0,  j ? 1 (4.20) 

The  square  of  the  ith  process  is 

Ni(t)  = Pi(t)  Gi(t)  (4.21) 

and  its  average  value  is 

E[N*(t)]  - (4.22; 

Higher  order  moments  are  evaluated  in  Section  5.0. 

4.3  Power  Envelope  Processes 

This  section  introduces  power  envelope  processes  that  can 
be  employed  in  the  class  of  compound  random  processes  discussed 
in  Section  4.2.  Statistical  properties  including  lower-order 
moment  functions  are  derived. 
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The  first  power  envelope  process,  which  may  be  termed  a 
normalized  chi-square  process,  is  an  obvious  extension  of  the 
chi-square  random  variable : 


P(t)  = pn-1  E G2(t) 
i=l  1 


(4.23) 


where  p>0  is  a constant  with  dimensions  of  power;  G^(t)  is  a 
stationary  Gaussian  process  with  zero  mean  and  unit  variance, 
statistically  independent  of  Gj(t),  j ^ i.  It  is  clear  from 
(4.23)  and  the  stipulation  above  that  P(t)  > 0.  Given  the 
properties  of  the  Gaussian  processes  G^(t)  stated  above,  the 
mean  value  of  P(t)  is 

nip  - p (4.24) 

The  second  moment  function  of  P(t)  is 

E[P(u)P(v)]  = (pn-1)2  E 2 E[G2(u)G2(v)]  (4.25) 

i-1  J-l  1 0 

For  i ? J,  the  fourth-order  moment  function  in  the  summand  is 

E[Gj(u)]EIG2(v)J  =1x1  (4.26) 

For  i = j,  the  moment  function  is  expanded  via  (7*28)  of  Ref.  12: 

E[G2(u)G2(v)]  = 1 + 2r2 (u-v)  (4.27) 

where  r(  ) is  the  autocovariance  coefficient  function  for  G1(t), 
all  i.  Substituting  (4.26)  and  (4.27)in  (4.25)  yields 
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E [P(u )P( v ) ] 


(pn-1  )2 


n n n 

£ £ 1 + £ [1  + 2r2(u-v)] 

i-1  J-l  i=l 


(4.28  ) 


Consolidating  the  first  term  in  the  single  summand  with  the  double 
sum  gives 


E[P(u)P(v)]-  - (pn-1)2  [n2  + 2nr2(u-v)]  (4.29) 

Then  subtracting  the  square  of  the  mean  gives  the  autocovariance 
function  of  P(t ) : 

Kp(u-v)  = 2p2n“ lr 2 (u-v)  (4.30) 

Evaluating  (4.3 0)  for  u = v and  dividing  by  mp  gives  the  coefficient 
of  variation: 


cp  ■ op  t mp  ■ ^2  r n (4.31) 

The  sequence  of  values  that  can  be  assumed  by  cp  is  ^2,  1,  ^ 2/3, 
^1/2 , >|2/5 , ...  . For  small  values  of  n,  the  gaps  between  succes- 

sive values  of  cp  are  fairly  large;  therefore,  it  might  not  be 
possible  to  achieve  an  accurate  approximation  of  a desired  value 
of  cp  by  selecting  the  best  value  for  n. 

The  third  moment  function  for  the  power  envelope  process  is 

m3(u,v,w)  - E[P(u)P(v)P(w)] 

- (P  * n)J  £ £ £ E[G?(u)G?(v)g£(w)J  (4.32) 

i-1  J-l  k-1  i j k 
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The  triple  sum  is  decomposed  according  to  the  identity  of  the 
indices,  and  the  expected  values  are  evaluated  via  (7-28)  of 
Ref.  (12). 

!n  n n n n 

Z Z Zl+Z  Z[l+  2rz(u-w)] 
i=l  J=1  k=l  i-1  k=l 


i^J  Jrtc  &1 


i?4c 


n n n n 

+ Z Z [1  + 2r2(u-v)]  + Z Z [1  + 2r2(v-w)] 


i-1  J=1 

1*J 


J-l  k=l 
j^k 


n 


+ Z [1  + 2r2(u-v)  + 2r2(v-w)  + 2r2(w-u)  + 8r(u-v)r(v-w)r(w- 
1=1 


-u)j| 


(4.331 


Consolidating  the  sums  with  the  Identical  summands  gives 

!n  n 

n3  + 2 Z Z [r2(u-v)  + r2(v-w)  + r2(w-u)] 
i-1  J-l 


n ) 

+ 8 Z r(u-v)r(v-w)r(w-u)  > 

1*1  ) 


(4.34) 


The  third  Joint  central  moment  is 


Pyk(u,v,w)  = EttP.Cu)  - n^]  |P  (v)  - mj]  [Pk(w)  - m^]} 
• myk(u,v,w)  - KJk(v-w)  - Kjj(>mi) 

- \ Vu'v)  - Tft 
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Substituting  (4,24  ),  (4.30  ) and  (4.34)  in  (4.35)  yields 

y(u,v,w)  = (2p) 3n-2r (u-v)r ( v-w)r(w-u)  (4.36) 

3 

In  some  cases,  estimates  of  the  mean  and  autocovariance 
of  the  power  envelope  are  available.  The  former  can  be  used 
to  obtain  a value  for  p via  (4.24),  and  the  latter  can  be 
evaluated  at  u * v to  obtain  a value  for  n via  (4.30).  The 
autocovariance  coefficient  function  can  then  be  determined 
from  (4.30).  These  elements  determine  the  third  central 
moment  function  for  P(t)  per  (4.36). 

The  second  power  envelope  process  is  also  synthesized 
from  squares  of  independent  Gaussian  processes.  In  this  case, 
the  weights  in  each  term  are  not  equal,  and  there  are  an  infinite 
number  of  terms.  The  objective  is  to  derive  a non-negative  pro- 
cess for  which  both  mean  and  autocovariance  can  be  specified. 

This  objective  is  (generally)  not  met  by  the  chi-square  process, 
discussed  previously. 


The  process  is  represented  by 

P(t)  - p(l  - a)  E a1  G,2(t)  (4.37) 

i-0  1 

where  p > 0 is  a constant  with  dimensions  of  power 
0 < a < 1 is  a constant 

G1(t)  is  a stationary  Gaussian  process  with  zero-mean 
and  unit  variance,  statistically  independent  of 
Gj(t),  J / i. 


It  is  clear  from  (4.37)  and  the  stipulations  above  that  P(t)  > 0. 
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Given  the  properties  of  the  Gaussian  processes  G^Ct)  stated  above 
the  mean  value  of  P(t)  is 


mp  = p(l  - a)  E a 
r 1=0 


Since  the  indicated  sum*  equals  (1  - a)-1,  the  result  is 


mp  - p 


The  second  moment  function  of  P(t)  is 


(4.38) 


(4.39) 


00  00 

E[P(u)P(v)]  = p2  (l-a) 2 E E ai+J  E [G2 ( u)G2 ( v ) ] (4.40) 

i = 0 J =0  1 J 

For  i ? J,  the  fourth-order  moment  function  In  the  summand  is 


E[G2(u)]  E [Gj (v) ] -1*1 


(4.41) 


For  i * J,  the  moment  function  is  expanded  via  (7*28)  of  Ref.  (12) 


E(G2(u)G2(v)]  = 1 + 2r2(u-v) 


(4.42) 


where  r(u-v)  is  the  autocovariance  coefficient  function  for  G^(t) 
all  i.  Substituting  (4.42  ) and  (4.41  ) in  (4.40  ) yields 


•Ref.  13,  0.231,  p.  7. 
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E[P(u)P(v)]  = p2 (1-a2 ) 
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ai+J  + E a2i  [1  + 2r2 (u-v) ] 
1=0 


(4.43) 


Consolidating  the  first  term  in  the  single  summand  with  the  double 
sum  yields 


E[P(u)P(v)]  = p2 (1-a2 ) 


E a1^  + 2r2 (u-v)  E a21j 


(4.44  ) 


Then  evaluating  the  sums  via  0.231  No.  1 of  Ref.  13  and  cancelling 
common  factors  gives 


E [P(u)P(v) ] = p2  [l  + 2 iz|r2(u-v)] 


(4.45) 


Utilizing  (4.45)  and  (4.39)  gives  the  autocovariance  function  of 
P(t)  as 


Kp(u-v)  = 2p2  i=|  r2 (u-v ) 


(4.46) 


Evaluating  the  function  for  u = v gives  the  variance  of  P(t)  as 


o|  = 2p2  (l-a)/(l+a) 


(4.47) 


Thus 


Kp(u-v)  = [opr(u-v)] 2 


(4.48) 


Using  (4.47)  and  (4.39)  gives  the  coefficient  of  variation  as 


cp  = op  * mp  = ^2(l-a)/l+a) 


(4.49) 
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It  is  seen  that 


lim  cp  = J2  (4 

a+0  ^ 


and 


lim  cp  = 0 (4 

a+1 


By  selecting  an  appropriate  value  for  the  constant  a,  the  coeffi- 
cient of  variation  of  P(t)  can  be  any  value  between  zero  and  iF 
For  a = 1,  the  sum  has  but  one  non-zero  term;  i.e.,  it  is  a chi- 
square  process  with  one  degree  of  freedom.  And  as  a + 1,  the 
process  becomes  deterministic.  Solving  (4.49)  for  the  constant  a 
yields 


a 


2 + c; 


(4.52) 


Given  the  coefficient  of  variation  p as  a process  parameter,  the 
associated  value  of  the  parameter  a can  be  determined.  For  the 
Kurtosis  of  the  compound  random  process  with  the  power  envelope 
process  given  by  (4.37),  substituting  (4.50)  and  (4.51)  in  (4.n) 
s -.-t 

1 lim  ya  t a"  - 9 (4.53) 

i a-0  4 P 


lim  Mb  * “ 3 (4.54) 

a-*-l  4 r 
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For  a-KL,  the  compound  process  degenerates  to  a stationary 
Gaussian  process. 

The  third  moment  function  for  the  power  envelope  process 
is 


m3(u,v,w)  = E [P(u)P(v)P(w) ] 


p (1-a)  E E i a 
i-0  j-0  k=0 


I ai+J+k 


E[Gi(u)Gj(v)Gk(w)]  (4*55) 


The  triple  sum  is  decomposed  according  to  the  identity  of  the  indices, 
and  the  expected  values  are  evaluated  via  (7.28)  of  Ref.  (12). 


m3(u,v,w) 


» p 3 (1-a) 3 < E 

00 

E 

E ai+J+k  + 

00  I 

E 

( i-0 

J-o 

k=0 

i-0  k 

i*J 

J^k 

k^i 

i/k 

OD  00  « . 

+ E E ai+2J 

[1  + 

2r2 (u-v) ] + 

00  00 

E E 

2i+k 


[1+  2r2  (u-w)  ] 


i-G  J-0 
i*J 


j-0  k-0 
J^k 


+ E a31  [1  + 2r2(u-v)  + 2r2(v-w)  + 2r2(w-u)  + 8r(u-v)r(v-w)r(w-u)]  | 
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Reconsolidating  the  sums  with  identical  summands  gives 


m_(u,v,w)  = p * (1-a)3  \ I Z a 

I \ 1*0 


OO  OO 

+ 2 Z a Z a2^  [r2(u-v)  + r2(v-w)  + r2(w-u)] 
1*0  J-0 


+ 8 Z r(u-v)r (v-w)r(w-u) 
i-0 


(4.57) 


Evaluating  the  indicated  sums  via  0.231  No.  1 of  Ref.  (13)  and 
cancelling  common  factors  yields 


m3(u,v,w)  * p3  | 1 + 2j=|  lr2 (u-v)  + r2(v-w)  + r2(w-u)] 

+ r(u-v)r (v-w)r(w-u)  j 

The  third  Joint  central  moment  is 
yijk(u,v’w)  5 EUP1(u)  - mj  IP  (v)  - mj]  [Pk(w)  - mk]} 

■ m3(u,v,w)  - m^Kj^v-w)  - mjKiic(w“u) 

- \Kij  <u-v>  - 


<1.58  ) 


(4.59) 


* 
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Substituting  (4.58  ),  (4.46  ),  and  (4.39)  in  (4.59)  yields 

y(u,v,w)  = $p3r(u-v)r(v-w)r(w-u ) (4.60  ) 

where  B = [2(l-a)]3  (1-a3)-1 


In  some  cases,  there  are  experimental  data  concerning  the 
mean  value  and  the  autocovariance  function  for  short-term  average 
power,  which  approximates  P(t).  These  data  can  be  used  to  es- 
tablish a value  for  p via  (4.39),  a value  for  op  and  the  auto- 
covariance function  r(  ) via  (4.48).  Prom  p and  ap  the  value  of 
the  weighting  parameter  a is  determined  via  (4.52).  These  elements 
determine  the  third  central  moment  function  per  ( 4 .6 0 ) . 


The  process  defined  by  (4.37)  meets  the  objectives  stated  in 
the  first  paragraph  of  this  section,  with  the  limitation  that  the 
variance  cannot  be  greater  than  twice  the  mean  value. 


Two  power  envelope  processes  have  been  proposed  and  analyzed. 
Given  that  parameter  values  are  selected  so  that  both  have  the 
same  mean  and  variance,  how  are  their  thiTd  central  moments  related? 
If  means  and  variances  are  equal,  the  coefficients  of  variation 
are  equal.  Equating  expressions  for  the  latter  and  solving  gives 
a = (n-1)  t (n+1).  Substituting  that  result  in  (4.60)  for  u * v 
gives  the  third  central  moment  for  the  WIS  (weighted  infinite  sum) 
process  in  terms  of  the  NCS  (normalized  chi-square)  parameter  n. 
Dividing  that  result  into  (4.36  ) evaluated  at  y * v and  simplify- 
ing gives 


u3(ncs)  . x 

U3(WIS)  “ t + TiHT 


(4.61) 
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Except  for  the  case  n = 1,  when  the  processes  are  identical,  the 
WIS  density  is  more  skewed  than  that  of  the  NCS. 

A set  of  power  envelope  processes,  some  or  all  of  which  may 
be  mutually  dependent,  can  be  synthesized  from  a set  of  independent 
power  envelope  processes: 


pl<t)  - | aa  yt) 


(4.62) 


where  a^  >_  0 


Q^(t)  > 0,  and  Q^t)  is  statistically  independent  of  Q.(t) 

mrt.  1 


The  mean  value  of  P^(t)  is 


where  m^  * E[Q^(t)] 


Z a^j 


L 


(4.63) 


The  Joint  second  moment  function  is 


Rij(u-v)  - E CPiCt )Pj (t ) 1 


Substituting  (4.62)  in  (4.64)  and  evaluating  gives 


R1J(U'V)  ' * * alm  ajn  EIVU>«„<T)1 

m n 


(4.64  ) 


(Z  a._  m )a  + Z a.  a,_  K_  (u-v)  (4.65) 

m im  m m lm  Jm  Qm 
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where  Kp.  ( ) is  the  autocovariance  function  for  Q (t).  Utiliz- 
Wm  m 

ing  (4.63)  and  (4.65)  gives  the  cross-covariance  function  for 
P±(u)  and  Pj(v): 

Kij(u-V)  ■ l aim  ajm  KQm(u-v)  <“•  66> 

It  is  seen  that  the  processes  P^(t)  and  Pj  ( t ) will  be  dependent 

if  a.  >0  and  a.  >0  for  at  least  one  value  of  the  index  m. 
im  j m 

The  third  joint  central  moment  function  is  evaluated  by 
means  of  a similar  procedure,  and  the  result  is 

Mijk(u>v,w)  = z alm  aJm  akm  yQra(u,v,w)  (4.67) 

where  yQrji(u>v»w)  is  the  third  central  moment  function  for  the 
power  envelope  process  Qra(t). 

If  the  basis  processes  are  either  of  the  types  discussed 
previously  in  this  section,  then  the  cross-covariance  function 
can  be  evaluated  by  substituting  (4.48)  in  (4.66): 

KiJ<u-v)  = Z alm  [%i-m(u-v)]»  (,.68) 


A similar  procedure  is  followed  to  obtain  the  third  joint  central 
moment  function. 
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ViiJk(u,v>w)  = Z a 

TP 


im  aJm  akm  ®mPmrm^u”v^rm(v-w)^ni(w“u)  (4-69  ) 


where  3m  = 2 3 n~z  for  the  envelope  process  defined  by  (4.23) 


[2(l-am)] 3 t (1-a^)  for  the  envelope  process  defined 


by  (4.37) 
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5.0  EVALUATION  OF  OUTPUT  MOMENTS 

Section  3*0  established  the  need  for  calculating  the 
moments  of  the  test  statistic.  The  moments  of  the  test 
statistic  to  order  three  will  be  evaluated  for  the  case  in 
which  the  inputs  to  the  multichannel  analog  are  compound 
random  processes  of  the  type  described  in  Section  4.0. 

The  moments  are  obtained  by  deriving  the  expected 
values  of  the  appropriate  powers  of  the  test  statistic  as 
given  by  (2.9),  which  can  be  expressed  in  more  compact  from 
as 


T f n 

Z * T-1  / du  2No(u)Nn(u)  - Z c 
0 L s 0 i=S 

where  Ng(t)  = S(t),  c^  = -1,  and  the  indices  in  the  sum  are 
S,  0,  1,  2,  3, • • • , n. 

The  first  moment  or  mean  value  of  the  test  statistic  is 
the  expected  value  of  (5.1): 


±N*(u) 


(5.1) 


- T-1  / du  | 2e[ns(u)Nq(U) ] - E ciE[Ni(u) ] } (5-2) 


The  first  term  of  the  integrand  is  zero  per  (4.19),  and  the 
components  of  the  indicated  sum  are  evaluated  per  (4.22). 
Evaluation  of  the  integral  then  yields 


n 

m7  » - Z c.p. 
Z i«S  1 1 


(5.3) 
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The  second  moment  is  the  expected  value  of  the  square 
of  (5-1),  which  can  be  expressed  as: 


T T 

Z2  = T"2  / du  / dv 
0 0 


j^4Ns(u)Ns(v)N0(u)N0(v; 


- 4Ns(u)Nq(u)  Z ( v)  + Z Z ciCjN^(u)Nj (v)j 


(5.4) 


If  the  first  term  of  the  integrand  is  expressed  in  terms  of 
(4.18),  its  expected  value  is 


E(TX)  -4E^Ps(u)Ps(v)P0(u)P0(v)]  e[gs(u)Gs(v)]  e[g0(u)G0(v)J  (5-5) 

and  if  the  relaxation  time  of  the  P-processes  is  much  longer 
than  those  of  the  G-processes,  then 


E(T1)  * 4psp0ps(u-v)p0(u-v) 


(5.6) 


If  the  second  term  of  the  integrand  is  expressed  in  terms  of 
(4.18),  its  expected  value  can  be  expressed  as 

E(T2)  - -4  JsE[^PS(u)P0(u)Pi(v)]  E[Gs(u)G0(u)G^(v)] 


(5.7) 


For  i"S,  the  second  expected  value  i3  E[Gs(u)Gg(v)]  E[GQ(u)l, 
which  is  zero  since  the  second  factor  is  zero.  For  i»0,  the 
second  expected  value  is  E[Gs(u)]  E [GQ(u)G^(v) ] , which  is  zero 
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' 


because  the  first  factor  is  zero.  And  for  i = 1,  2,  n, 

the  second  expected  value  is  E[Gs(u)]  E[Go(u)]  E[G2(v)],  which 
also  equals  zero.  Thus 


E(T2)  * 0 (5.8) 

The  third  term  of  the  Integrand  of  (5.4)  is  expressed  as  a 
double  summation.  For  i ^ J, 

E[N2(u)N2(v)]  = E[P1(u)PJ(v)]  E[Q2(u)]  E[Q2(v)] 

* PjPj  + Kij(u-v)  (5.9) 

since  both  the  second  and  third  factors  are  unity.  For  terms 
such  thaw  i ■ J 

E[N2(u)N2(v)]  - E[P1(u)P1(v)]  E[Q2(u)Q2(v)]  (5-10) 

The  second  factor,  a fourth  moment  function  of  a Gaussian  process, 
is  evaluated  by  means  of  (7.28)  of  Ref.  12.  Thus 

E[Nj(u)N*(v)J  - [p2  + K11(u-v)]  [1  + 2p2 (u-v) ] (5.11) 

“ pi  + Kii^u“v)  + 2 (p^  + o2)p2(u-v)  (5.12) 

Utilizing  (5*9)  and  (5-12)  permits  the  evaluation  of  the  expected 
value  of  the  third  term  of  the  integrand: 

n n n 

E(T,)  - I Z c.c,  [p.p.  + K. , (u-v) ] + 2 Z c?(p?  + o2 ) p2(u-v) 
j i*S  J*S  1 J 1 J i"S  x A 1 

(5.13) 
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where  W"1  = 2 / dyp,(y)p,(y)  can  be  considered  as  the 
1J  _co  1 J 

reciprocal  of  the  joint  equivalent  bandwidth 
of  the  associated  random  processes. 


The  third  moment  of  the  test  statistic  is  the  expected 
value  of  the  cube  of  the  test  statistic,  which  can  be  expressed 


T T T r 

Z 3 = T“3  / du  / du  / dw  8Nq(u)N-(v)NQ(w)Nn(u)Nn(v)Nn(w; 

oooLSbSOO° 


- 12Ns(u)Ns(v)N0(u)N0(v)  Z C^^w) 

i=S 


n n 

+ 6N„(u)Nn(u)  Z Z c,c,N^(v)N3(w) 
S 0 i-S  J-S  1 J 1 J 


n n n -| 

- Z Z Z c.c  .c.N.2(u)N2(v)N£(w) 

i-S  j-S  k=S  1 J K 1 J K J 


(5.18) 


If  the  first  term  of  the  integrand  is  expressed  in  terms 
of  (4.18),  its  expected  value  is 

E(UX)  - 8E  ^Ps(u)Ps(v)Ps(w)P0(u)P0(v)P0(w)] 

x e[qs(u)Gs(v)Gs(w)]  E [g0(u)G0(v)G0(w)]  (5.19) 

According  to  (7-28)  of  Ref.  12,  the  second  and  third  factors  are 
zero;  thus 


E(UX)  - 0 


(5.20) 
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If  the  second  term  of  the  integrand  of  (5-18)  is  expressed 
in  terms  of  (4.18)  and  then  partially  expanded,  its  expected 
value  is 

E(U2)  = 12 |e  ^Ps(u)Ps(v)P|(w)P0(u)P0(v) 
x E |gs(u)Gs(v)G|(w)J  E^G0(u)G0(v)J 

+ E [ ^Ps(u)Ps(v)P0(u)P0(v)p£(w)  ] 
x E^Gs(u)Gs(v)J  E ^Gq  (u)Gq ( v )Gq(w)J 

- c1e[i|pS(u)PS(v)P0(u)P0(v)P*(w)] 

x e[gs(u)Gs(v)]  e[g0(u)Gq(v)]  e[g*(w)]|  (5.21) 

Applying  (7.28)  of  Ref.t  12  yields 

E(U2)  * 12  |e  ^Ps(u)Ps(v)P0(u)P0(v)P|(w) 


x p0(u-v)|^ps(u-v)  + 2ps(w-u)ps(v-w) 


+ E [^Ps(u)Ps(v)P0(u)P0(v)P5(w)] 
x Ps(u-v)  j^p0(u-v)  + 2p0(w-u)p0(v-w)J 

" 1f1  ciE  [l|PS(u)PS(v)P0(u)P0(v)Pi(w)] 
x ps(u-v)p0(u-v)| 


(5.22) 
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Combining  the  first,  third,  and  fifth  terms  gives 


E(U  ) = 12 
2 


| 2E  ^Ps(u)Ps(v)P0(u)P0(v)P|(w)  ] 


pQ(u-v)ps(w-u)ps(v-w) 


2E  ^Ps(u)Ps(v)P0(u)P0(v)P*(w)  j 


x ps(u-v)p0(w-u)p0(v-w) 


- Js  c±E  [^Ps(u)Ps(v)P0(u)P0(v)Pj(w)] 
x ps(u-v)p0(u-v)  | 


(5.23) 


If  the  relaxation  times  of  the  P-processes  are  much  longer  than 
those  of  the  G-processes, 

E(U2)  « 12  j 2E^P|(w)P0(w)Jps(v-w)ps(w-u)p0(u-v) 

+ 2E jps(w)P*(w)J  ps(u-v)p0(v-w)p0(w-u) 


- ^ c1e|ps(u)P0(u)P1(w)^ps(u-v)p0(u-v)| 


(5.24) 


If  the  third  term  of  the  integrand  of  (5.18)  is  expressed  in  terms 
of  (4.18),  its  expected  value  is 


Bolt  Beranek  and  Newman  Inc. 


Report  No.  3656 


E(V  " 6 ^3  j£g  ci°jE[^PS(u)P0<u)Pi(v>PJ(w)] 

X E jGs(u)GQ(u)G|(v)Gj(w)J 


(5-25) 


For  terms  In  which  neither  i nor  j is  equal  to  S or  0,  the  second 
expected  value  is  E [Gs(u) ]E [G0(u) ]E [Gi(v)Gj (w) ] = 0.  For  i = j = s, 
and  for  i = J * 0,  the  second  expected  value  in  (5-25)  is  also  zero. 
Finally,  for  i = S and  J = 0 and  vice  versa,  the  second  expected 
value  Is  also  zero.  Thus 


E(U3)  = 0. 


(5-26) 


The  fourth  term  of  the  integrand  of  (5-18)  is  represented  as 
a triple  sum.  Within  this  sum  there  are  n + 2 terms  with  i = J = k. 
The  expected  value  of  the  sum  of  those  terms  can  be  represented  as 

E(U41)  = - Z c^E^P1(u)P1(v)P1(w)Je[g^(u)G^(v)G^(w)]  (5.27) 

Applying  (7.28)  of  Ref.  12  to  the  second  expected  value  gives 

E(UU)  * - ciE[p1(u)P1(v)P1(w)J  £ 1 + 2p^(u-v) 

+ 2p*(v-w)  + 2p^(w-u)  + 8p1(u-v)p1(v-w)pi(w-u)J  (5.28) 

In  the  triple  sum  there  are  3(n+2)(n+l)  terms  in  which  two  of  the 
Indices  are  equal  but  different  from  the  third.  After  application 
of  (7.28)  of  Ref.  12,  the  expected  value  of  the  sum  of  those  terms 
can  be  represented  as 
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E(U«2>  * - Js  sls  alCj{E[PlCu)PlCv)PJ  CW)]  t1  + 


+ E 
+ E 


Pj(u)P1(v)P1(w) 

pi ( U. ) Pj  ( V 


1 + 2p^( v-w) 


1 + 2pi(w-u) 


(5-29) 


In  the  triple  sum  there  are  (n+2)(n+l)n  terms  with  all  of  the 
summation  indices  different.  The  expected  value  of  the  sum  of 
these  terms  can  be  represented  as 


n n n r 1 

E(UaO  = - E E Z c.c.c.Elp  (u)P.(v)P,(w)I 
i=S  j=S  k=S  1 J K L 1 J J 


(5.30) 


J 

i^J  j^k  kj^i 


Summing  (5.28)  through  (5-30)  gives  the  expected  value  of  the 
fourth  term  of  the  integrand  of  (5.18): 

E(Ua)  = - E E E c.c.c  e[p.(u)P.(v)P,  (w) 

4 i=S  j=S  k*S  1 J K L 1 J J 

■ 2 is  Js  =2i<:j{E[pi(u)pi<v>pj<',>]‘jii<u-v) 

+ E jjj  (u)P1(v)P1(w)J  p^(v-w)  + eJpi(u)Pj  (v)Pi(w)j  p^w-u)! 
:31E|pi(u)P1(v)P1(w)Jp1(u-v)p1(y-w)  Pl(w-u) 


n 

- 8 E 


(5.3D 
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Expressions  (5-24)  and  (5.31)  give  the  expected  values  of 
the  terms  of  the  integrand  of  (5.18).  Given  the  symmetry  of 
certain  of  the  integrands,  the  expected  value  of  the  cube  of 
the  test  statistic  can  be  expressed  as 

T t T f 

E(Z*)  = T“*  / du  / dv  / dw  24E["Po(w)Pr.(w)l  pQ(v-w)p  (w-u)pn(u-v) 
000  L&UJ^  S u 

+ 24E^Ps(w)P*(w)J  Ps(u-v)pQ(v-w)p0(w-u) 

- 12  ^E^  c1E^Ps(u)P0(u)P1(w)j  ps(u-v)pQ(u-v) 

n r 1 

- 8 z c^E|^P1(u)P1(v)P1(w)J  p1(u-v)p1(v-w)p1(w-u) 

- 6 ^E^  ^ E^  c^CjE  jV1(u)P1(v)Pj  (w)J  p^(u-v) 

n n n r . ) 

- E E E c.c.c  E P (u)P.(v)P.  (w)  (5.32) 

i-S  J-S  k-S  1 J K L 1 J K J ) 

The  third  central  moment  of  the  test  statistic  is 


- E(Z-mz)» 


E(Z»)  - 3m zo*  - mj 


(5.33) 


The  third  moment  functions  appearing  in  the  integrand  can  be 
expressed  in  terms  of  central  moment  functions.  The  most 
general  of  these  is 
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E [pl(u)Pj  (v)pfc(w)] 

■ E{[pi<u)  - pl + Pi]  [pJ(v)  - Pj  + Pj]  [pk(w) 

■ y1Jk(u-v,v-w,w-u)  + PiKjjjCv-w)  + PjKki(w-u) 

♦ PkKu(u-v)  + Pipjpk 


pk + pk. 


(5.34) 


where  p1Jk(u-v,v-w,w-u)  = EUP±(u)  - p±]  IPj(u)  - p ^ ] [Pk(u)  - pkJ}, 

The  expression  for  E(Z3)  is  expanded  by  evaluating  (5*34)  for  the 
various  terms  of  (5*32).  The  expanded  result  is  substituted  in 
(5*33),  along  with  (5-3)  and  (5.14)  to  obtain  the  third  central 
moment  of  the  test  statistic,  shown  in  Table  1. 


Table  1.  Third  Central  Moment  of  Test  Statistic 


y3z = r>  \ du  l dv  \ dw  1 24[psso  + 2Ps°so  + Po(ps+as)]ps{v-w)ps(w-u,po(u-v) 

+ 24[psoo  + 2poaso  + ps(po+ao)Jps(u'v)po(v-w,po(w-UJ 

- 12  Z$  cijySQ1(0,u-w,w-u)  + psK01(u-w)  + PqKj^w-u)  + P^^jpjtu-vJpQtu-v) 


- 8 cj  ju31  + Sp^oJ  + pjjp^u-vjp^v-wjp^w-u) 

- 6 Z$  Z$  cjcj ^^(O.u-w.w-u)  + 2piK1j(u-w)jpj(u-v) 


n n n 


Z Z Z CjC.c. y. ik(u-v,v-w,w-u) 
1-S  J-S  k-S  1 J k 1JK 
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The  first,  second,  and  fourth  lines  in  Table  1 are  repre- 
sented by  an  integral  of  the  form 


T/2  T/2  T/2 

I * T “3  / du  / d v / dw  r(u-v)s(v-w)t(w-u)  (5.35) 

-T/2  -T/2  J-T/2 


A change  of  variables  x * w-u  and  y ■ w+u  yields,  for  half  of  the 
domain  of  integration 


T/2  T T-x 

I + “ \ T“ 3 / dv  / dx  t(x)  / dy  r(y/2-x/2-v)s(v-y/2-x/2)  (5-36) 

J-T/2  0 -T+x 


which  equals  the  result  from  the  other  half  of  the  domain  of  inte- 
gration, given  that  all  of  the  functions  have  even  symmetry.  If 
the  decay  time  of  t(x)  is  small  compared  to  the  integration  period 
T,  then 


T/2  T T 

1*21+  * T-3  / dv  / dx  t(x)  / dy  r(y/2-x/2-v)s(v-y/2-x/2)  (5-37) 

-T/2  0 m 


-T 


If  the  functions  r(  ) and  s(  ) are  represented  as  the  inverse 
Fourier  transforms  of  R(  ) and  S(  ) respectively,  then,  after 
some  algebraic  manipulation. 


I - T-»  f df  R(f)  f dg  S(g)  f dv  exp  2irj(g-f)v 
-•  -oa  -T/2 


Tr  T 

xj  dx  t(x)  exp  I-irJ  (f+g)x]  / dy  exp  irj(f-g)y 

A 


(5.38) 
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* 

\ ' 


If  the  decay  time  of  t(  ) is  small  compared  to  T,  then  the 
fourth  integral  of  (5-38)  is  approximately  where  T(  ) 

is  the  Fourier  transform  of  t(  ).  Substituting  this  result  in 
(5-38)  and  evaluating  the  v and  y integrals  yields 


= T"1  f 


L 


df  R(f)  / dg  S(g)  T 


. 2 

/f+gl  /sin  ir(g-f )t\ 
V 2 / (g-f)T  ) 


(5.39) 


A change  of  variables  h = g - f in  the  second  integral  yields 


00 

* / 

2 

\ / \ 

/ df  R(f) 

im  00 

f dh  S(h+f ) TMj+1 

fsin  irhTi 

J y uhT  ) 

(5.40) 

If  the  averaging  time  T is  large,  then  the  last  factor  is  very 
small  except  in  a small  increment  around  h = 0.  If  S(  ) and  T(  ) 
are  bandpass  functions,  then  f>>h  over  the  bandpasses,  and 


I 


T"1 


00 


df  R(f ) S(f)  T(f ) 


(5.41) 


The  h-integral,  evaluated  via 3- 821,  number  9 of  Ref.  13,  equals 
T_1j  thus 


I = T~ 


df  R(f ) S(f)  T(f ) 


(5.42) 


Since  the  spectral  functions  in  (5.42)  are  Fourier  transforms  of 
autocovariance  coefficients,  the  dimension  of  the  integral  in 
(5.42)  is  1 * frequency  squared.  Let 


I 
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Krt) 


(5.44) 


The  third  and  fifth  lines  of  Table  1 include  integrals  of 
the  form 


where  the  decay  times  of  r(  ) and  s(  ) are  much  smaller  than  that 
of  f(  ),  and  much  smaller  than  the  length  of  the  integration 
period  T.  A change  of  variables  x * u-v  and  y ■ u+v  yields,  for 
half  of  the  domain  of  integration. 


I + 


T T/2  T/2 

h T ”3J  dx  r(x)s(x)  f dw  f dy  f(x/2+y/2-w)  (5.46) 

* : -T/2  -T+x 


/ 


Wrft  = 4 I df  R(f)  S(f)  T(f) 


(5.43) 


If  all  of  the  spectra  have  coincident  rectangular  bandpasses  of 
width  W,  then  Wrft  » W.  Utilizing  (5.43)  in  (5.42)  gives 


T/2  T/2  T/2 

I * T~  / du  / dv  / dw  f (u-w)r(u-v)s(u-v) 
-T/2  -T/2  ^T/2 


3 


(5.45) 


which  equals  the  result  from  the  other  half  of  the  domain  of  in- 
tegration. Given  the  stipulations  stated  above 


I 


I 

I 
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A change  of  variables  to  p * y/2-w  and  q = y/2+w  gives,  for 
half  of  the  domain  of  integration 


I » T-3  / dx  r(x)s(x)  / dp  f(p)  / dq, 
v 0 0 -T+p 


which  equals  the  result  from  the  other  half  of  the  domain. 
Evaluating  the  q integral  yields 


I ■ 2T  * 4T“2  / dx  r(x)s(x)  f dp  f(p)  (1-pT-1) 


Given  the  stipulations  stated  above. 


I « 2T“2  / dx  r(x)s(x)  f dp  f(p)  (l-pT“l) 

-OO  0 


(2Tl«rS)" 


/ dp  f(p)  (1- | p | T“l ) 
-T 


(5.48) 


(5-49) 


(5.50) 


(5-51) 


where  the  Joint  equivalent  bandwidth  is  defined  below  (5.17). 

Application  of  (5-46)  and  (5-51)  to  the  expression  in  Table  1 
yields  the  result  in  Table  2. 

The  last  term  of  the  expression  in  Table  2 includes  a triple 
Integral  whose  Integrand  is  the  Joint  third  central  moment  func- 
tion of  the  power  envelope  processes.  This  integral  can  be  partially 
evaluated  for  the  power  envelope  processes  defined  by  (4.23)  and 
(4.37).  For  those  processes,  the  Joint  third  central  moment  function 


I 

I 
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i 

) 

f 


is  given  by  (4.36)  and  (4.60)  respectively;  for  either  case, 
the  essential  integral  can  be  represented  by 

T T T 

I = / du  / dv  / dw  r(u-v)s(v-w)t(w-u)  (5.52) 

0 0 0 

A new  set  of  variables  is  defined  by  the  transformations 

u = x - z x = -v+w 

v»  y - z y = -u  +w 

w * x + y - z z = -u  -v  +w 


With  the  change  of  variables,  the  Integral  is 

I = / dx  / dy  / dz  r (x-y )r (-x)r (y ) 

R 

The  region  of  integration  is  a three-dimensional  polygon. 

Figure  5.1  shows  a plan  view  of  the  polygon  in  the  x-y  plane. 
The  z-coordinate  of  each  vertex  is  indicated.  The  figure  also 
gives  the  equations  of  the  faces  of  the  polygon.  The  region  of 
integration  is  divided  into  six  parts  numbered  as  shown. 

Part  1 of  the  region  is  bounded  above  by  the  face  ABCD  and 
below  by  the  face  CDHO.  It  is  also  bounded  by  the  plane  x*y, 
and  y-0.  The  integral  for  this  part  is  therefore 


(5.53) 


(5.54) 
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0 0 

f dx  / dy  r(x-y )r(-x)r(y) (T+x) 
-T  x 


(5-56) 


An  alternative  form  is  found  by  changing  the  signs  of  the  variables 
of  integration: 


T x 

I,  * / dx  / dy  r(-x+y )r(x)r(-y) (T-x) 
1 0 0 


(5-57) 


Part  2 of  the  region  is  bounded  above  by  the  face  ADHE  and 
below  by  DCHG.  It  is  also  bounded  by  the  planes  y=0  and  x=0. 
The  integral  for  this  part  is 


T -T+y  x 

I?  = / dy  / dx  r(x-y)r(-x)r(y)  / dz 
^00  -T+y 


T -T+y 

f dy  f dx  r (x-y)r(-x)r(y) (T+x-y ) 
0 0 


Now  let  w ■ -x+y  in  the  inner  integral: 

T T 

I5  = / dy  / dw  r (-w)r (w-y)r (y) (T-w) 

* 0 y 

Interchanging  the  order  of  integration  then  gives 

T w 

I-  * / dw  / dy  r(-w)r(w-y)r(y) (T-2) 

£ 0 0 


(5.58) 


(5.59) 


(5.60) 


(5.61) 


I 
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The  function  r(  ) is  the  autocovariance  coefficient  function  for 
the  Gaussian  processes  that  are  the  components  of  the  power  en- 
velope processes  described  in  Section  4.3;  thus,  it  is  an  even 
function  of  its  argument,  and  comparing  (5.61)  to  (5*57)  shows 
that  I2  = I . 

Part  3 of  the  region  is  bounded  above  by  the  face  ADHE  and 
below  by  EFGH.  It  is  also  bounded  by  the  planes  x = 0,  and  x = y. 
The  integral  for  this  part  is 


T y x 

/ dy  / dx  r(x-y)r(-x)r(y)  / dz 

0 0 -T+x+y 


T y 

/ dy  / dx  r (x-y )r (-x)r (y ) (T-y ) 
0 0 


(5-62) 


(5.63) 


Interchanging  the  variables  of  integration  in  (5.63)  and  comparing 
the  result  to  (5.57)  for  I,  shows  that  1^  = 1^. 

The  integrals  for  Parts  4,  5,  and  6 of  the  region  of  integra- 
tion are  identical  to  those  of  Parts  1,  2,  and  3,  respectively. 
Since  the  latter  trio  are  all  equal  the  complete  result  may  be 
expressed  by 


T x 

I - 6 / dx  r (x) (T-x)  / dy  r(x-y)r(y) 
0 0 


(5.64) 


Substituting  this  result  into  the  sixth  term  of  the  expression 
in  Table  2 gives 


wjjkZssmmt 
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n n n 

- Z 8 pi,  I Z Z Z c.a.  c.a.  c,  a. 
m m rm  m1=s  j_s  k=s  1 im  J Jm  k kra 


(5.65) 


where  1=6  T-2  / dx  r (x)  (1-xT-1)  / dy  r (x-y)r  (y) 


The  results  of  this  section  are  low-order  moments  of  the 
test  statistic  Z of  the  channel  of  interest.  The  results  include 
(5-3)  for  the  mean  value,  (5.17)  for  the  variance.  Table  2 for 
the  third  central  moment,  and  (5*65)  for  a partial  evaluation  of 
the  sixth  term  in  Table  2. 
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6.0  TIME-INVARIANT  POWER  ENVELOPES 


This  section  considers  two  cases  in  which  the  power 
envelopes  are  t ime -invariant . 


In  the  first  case,  the  envelopes  are  deterministic 
constants,  and  the  input  processes  are  therefore  Gaussian. 
This  case  is  analyzed  to  establish  a set  of  values  for  the 
weighting  coefficients  of  the  multichannel  analog. 


In  the  second  case,  the  envelopes  are  random  variables, 
and  the  input  processes  are  therefore  only  conditionally 
Gaussian  on  the  random  variables.  This  case  is  analyzed  to 
establish  the  conditions  for  the  applicability  of  a formula 
that  has  been  employed  to  calculate  the  probability  of  detec- 
tion in  an  environment  with  fluctuating  sonar  parameters. 

6.1  Special  Case:  Pj(t)  = p^ 

In  many  analyses  of  passive  sonar  receiver  performance, 
it  is  assumed  that  both  the  signal  and  the  noise  in  a single 
channel  are  zero-mean  statistically  independent  Gaussian  pro- 
cesses. For  the  analysis  of  this  section  it  will  further  be 
assumed  that  the  noise  processes  of  the  other  channels  are 
zero-mean  statistically  independent  Gaussian  processes,  all 
statistically  identical.  This  is  a degenerate  case  of  the 
compound  process  of  Section  4.0  in  which 


Pq(t)  = p< 


(6.1) 


P^($)  * Pjj*  1 * 0»  1>  2,  ...,  n 


(6.2) 
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I 

I 

I 


I 
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where  ps  and  pN  are  constants.  . The  results  of  this  case  will 
be  used  to  establish  a set  of  values  for  the  threshold  weight- 
ing coefficients  Cj^,  i = 1,  2,  3,  . ..,  n. 


The  mean  value  of  the  test  statistic  for  this  case  is 
found  by  substituting  (6.1)  and  (6.2)  in  (5*3): 


m 


Z 


(6.3) 


The  variance  of  the  test  statistic  is  found  by  substituting 
(6.1)  and  (6.2)  in  (5-17): 


2 


( WT ) “ 1 


PS  + 2pSPN  + PN 


n 

Z 

i-1 


(6.4) 


In  this  step,  it  was  assumed  that  the  equivalent  bandwidth  of 
the  signal  is  identical  to  those  of  the  noise  processes. 


The  probabilities  of  detection  and  false  alarm  are  given 
by  (2.13)  and  (2.14)  respectively,  which  show  that  both  of 
these  probabilities  depend  on  the  coefficient  of  variation 
CZ  5 9Z  * mZ  the  test  statistic.  The  inverse  of  that 
coefficient  for  the  case  of  signal  present  is  found  by  dividing 
(6.3)  by  the  square  root  of  (6.4): 


n 

r + 1 - Z c 
i-1 


♦ 2r  + 1 + 


n 

Z 

i-1 


(6.5) 


r 

* 

i 

i 

( 

i 


* 

i 


<a  'mwtWfr— r - 


Bolt  Beranek  and  Newman  Inc. 


Report  No.  3656 


where  r = pg  r is  the  signal-to-noise  power  ratio.  For 
the  case  of  noise  alone. 


n 


(6.6) 


This  result  is  independent  of  the  noise  power  pN- 

The  weighting  coefficients  will  be  selected  to  maximize 
the  detection  probability  while  achieving  a specified  false 
alarm  probability.  If  the  time-bandwidth  product  is  very 
large,  then  the  test  statistic  is  nearly  normal,  and  the  form 
of  the  density  function  for  the  standardized  variable  (Z  - mz) 
t oz  is  very  insensitive  to  the  choice  of  weighting  coefficients. 
Then  the  probabilities  of  detection  and  false  alarm  depend  al- 
most entirely  on  the  coefficient  of  variation  of  the  test 
statistic.  To  achieve  a specified  probability  of  false  alarm, 

l/c2  = -y  (6.7) 

where  y is  a positive  constant.  Substituting  (6.7)  in  (6.6) 
and  solving  gives 


and  substituting  this  result  in  (6.5)  gives 
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1_  m ^jWT  r - yq 

t _ 

^r2  + 2r  + q2 


(6.9) 


where 


q 


4 


n 

1 + E 
i=l 


(6.10) 


It  Is  clear  from  (2.13)  that  the  probability  of  detection  in- 
creases with  the  ratio  l/c^.  By  differentiating  (6.9)  with 
respect  to  q it  is  found  that  c^-1  decreases  monotonically  with 
q.  Thus  the  probability  of  detection  will  be  maximized  when 
q is  minimized.  Since  q is  invariant  with  an  interchange  of 
coefficients,  the  solutions  c^  that  minimize  it  all  have  the 
same  value,  say  c.  Substituting  that  symbol  in  (6.8)  yields 


^WT  (1  - nc ) = -y  ^1  + nc2  (6.11) 

Squaring  (6.11)  yields  a quadratic  equation  in  c,  and  substitut 
ing  the  appropriate  solution  of  that  equation  in  (6.10)  yields 


If 


a - 3 +lln(n  + 1 - Bz) 
n - B2 


(6.12) 


(6.13) 
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then 


q 


■4 


1 + n 


-l 


(6.14) 


It  Is  clear  from  (6.9)  that  the  signal-to-nolse  power  ratio 
required  to  achieve  a probability  of  detection  of  0.5  is 


rT  = Y Q * ^*WT 


(6.15) 


The  quantity  10  log10  rT  is  usually  known  as  the  recognition 
differential.  It  is  clear  from  (6.14)  that  performance  improves 
monotcnically  with  n.  However,  the  penalty  associated  with 
small  n is  not  severe;  for  example,  the  degradation  associated 
with  n = 2 is  only  0.9  dB. 


Figure  6.1  shows  transition  curves  for  cases  in  which 
the  bandwidth  is  600  Hz,  the  averaging  time  is  five  minutes, 
and  the  false  alarm  probability  per  channel  is  10-5.  The 
upper  plot  shows  detection  probability  versus  signal-to-noise 
power  ratio  for  the  case  n * 2 and  the  limiting  case  n ®, 
which  represents  perfect  estimation  of  background  noise  level. 
The  spread  of  the  curve  for  n ■ 2 is  noticeably  greater  than 


that  for  n -*■  « . (Spread  is  the  value  of  r for  Pn  * 0.8  less 


that  for  PD  = 0.2.)  To  a first  level  of  approximation,  spread 


is  proportional  to  q.  The  lower  plot  shows  detection  prob- 
ability for  the  same  cases  as  a function  of  10  logio  r.  To 
a first  level  of  approximation,  the  spread  of  these  curves  is 
independent  of  q. 


6.2  Special  Case: 


P^t) 


N 


This  section  extends  the  results  of  the  previous  section 
to  apply  to  cases  in  which  sonar  detection  experiments  are 
conducted  in  randomly  selected  environments.  The  analysis 
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is  based  on  a special  case  of  the  compound  process  in  which 

Ps(t)  = Ps  (6.16) 

pi(t)  - PN,  i - 0,  1,  2,  ...,  n (6.17) 

where  Ps  and  PN  are  statistically  independent  random  variables. 
The  signal  and  noise  processes  are  thereby  non-erogodic  and 
non-Gaussian;  however,  they  are  conditionally  Gaussian  given 
Pg  and  Pjj.  The  results  also  apply  to  cases  in  which  the  relaxa- 
tion times  of  the  power  envelope  processes  are  very  long 
compared  to  the  post-rectification  averaging  period. 

The  result  (6.9)  of  the  previous  section  is  utilized  to 
obtain  the  inverse  of  the  conditional  coefficient  of  variation 
of  the  test  statistic: 


, •Jwt  r - yq 

— = 1 (6.18) 

Z ^ R2  + 2R  + q2 

where  R ■ Pg  * Pn»  a random  variable,  is  the  signal-to-noise- 
power  ratio  for  a single  experiment. 

If  the  time -bandwidth  product  is  adequate,  the  conditional 
distribution  of  the  test  statistic  is  approximately  normal,  and 
the  conditional  probability  of  detection  can  be  expressed  as 

P(D|R)  - N( 1 ♦ Oj)  (6.19) 


r 

* 

r 

I 


i 


t 

( 
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where  N(  ) is  the  probability  distribution  function  for  a zero- 
mean  unit  variance  normal  variate.  The  dependence  of  0%  on  R 
is  given  by  (6.18).  The  probability  of  detection  is 

PD  = E[N(l/cz)l 


00 


* / dr  fR(r)  N 
0 R 


> r2  + 2r  + q2 


(6.20) 


where  fR(  ) is  the  probability  density  function  for  R,  which  can 
be  calculated  from* 


^(r)  = / ds  s fg(rs)fN(s)  (6.21) 

where  fg(  ) is  the  probability  density  function  for  Ps  and 
fN(  ) is  the  probability  density  function  for  PN. 

If  the  distribution  of  R is  broad  with  respect  to  the  width 
of  the  transition  curve,  then  the  latter  can  be  approximated  by 
a step  function  with  the  transition  occuring  at  rT,  defined  per 
(6.15),  and  the  approximation  of  (6.20)  is 

00 

PA  - / dr  f (r)  (6.22) 

rT 

A quantity  that  is  often  employed  in  sonar  analysis  is 

Ne  = io  log10  (R  * rT)  (6.23) 


*Ref . 5,  page  197. 
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which  is  the  so-called  "signal  excess",  the  signal-to-noise 
ratio  in  dB  in  excess  of  that  required  for  a conditional 
probability  of  detection  of  0.5.  If  the  corresponding  change 
of  variables  is  made  in  (6.22)',  the  result  is 


PA  « / de  fE(e)  (6.24) 

where  fg(e)  = 0. 23r<pe°*  ^3e  fR(rTe°*23e)  is  the  probability  density 
function  for  the  excess  signal-to-noise  ratio  in  dB.  The  form 
(6.24)  is  often  employed  with  the  assumption  that  E is  normally 
distributed.  Reference  14  compares  the  predictions  made  using 
that  form  with  NE  normally  distributed  to  those  using  the  density 
function  for  a short-noise  process. 


The  difference  between  the  predictions  of  (6.20)  and  its 
approximation  (6.22)  will  be  shown  by  means  of  several  examples. 
It  is  assumed  that  both  signal  and  noise  envelopes  are  gamma 
variables*  with  respective  densities 


fp(u) 


gp(u) 


p7W(^uj  e*p(-iu),u>0 

e,pUu)'ul° 


(6.25) 

(6.26) 


where  ps  - E(Pg) 


PN  - E(P„) 


*The  utilization  of  the  gamma  density  for  the  noise  envelope  is 
suggested  in  Reference  15. 
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The  parameters  a and  b determine  the  respective  coefficients  of 
variation  of  the  envelopes.  For  example,  the  coefficient  of 
variation  for  the  signal  is  1 * 

Substituting  (6.25)  and  (6.26)  in  (6.21)  and  evaluating 
the  integral  gives  the  density  of  the  signal-to-noise  power 
ratio : 


fR(r) 


_ f(a+b) 


ap 


f (a) T(b ) l bp{ 


N 


.a-1 


RFf 


(6.27) 


The  result  is  seen  to  be  a generalization  of  the  F density  in 
which  ps  t pN  is  a scale  factor.*  The  mode,  mean,  and  variance 
of  R are 


mode 


a-1 


b 

b+1 


N 


(6.28) 


1 

b PS 

mR  ■ b-1  * PN  * 

b > 1 

(6.29) 

’ 

. 1 

1 

1 

1 

_2  _ b+a-1  _2 

% sTb^iy  mR  ■ 

b > 2 

(6.30) 

III 

r 

•The  special  case  a * 1 corresponds  to  a Rayleigh  fading  (narrow- 
band)  signal.  This  case  was  discussed  in  a presentation  at  the 

94th  meeting  of  the  Acoustical  Society  of  America,  December,  1977, 
by  J.  C.  Heine  and  J.  R.  Nitsche. 

1 

. 

: 1 

. n 
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Dividing  (6.29)  by  the  square  root  of  (6. 30)  gives  the  coefficient 
of  variation 


CR  = 


b+a-1 
a(b-2 ) 


b > 2 


(6.31) 


It  is  seen  from  (6.29)  that  mR  > ps  * pN;  that  is,  ps  r pN  is 
not  the  mean  signal-to-noise  power  ratio. 

The  probability  of  detection  and  its  approximation  will  be 
calculated  by  substituting  (6.27)  in  (6.20)  and  (6.22)  respectively. 
The  receiver  parameters  are  those  of  one  of  the  cases  examined  in 
Section  6.1  in  which  the  bandwidth  is  600  Hz,  the  averaging  time 
is  five  minutes,  the  false  alarm  probability  per  channel  is  10-5, 
and  the  noise  background  is  estimated  with  two  channels.  Figures 
6.2  through  6.7  show  plots  of  PD  - PA  as  a function  of  Pp  for  six 
pairs  of  values  of  the  parameters  a and  b.  Peak  errors  range  from 
nearly  0.01  to  0.03.  For  detection  probabilities  greater  than 
about  0.2,  these  errors  are  small.  For  smaller  probabilities, 
the  percentage  of  error  increases. 

For  integer  values  of  the  parameters  a and  b in  the  density 
function  (6.27),  the  integral  (6.22)  can  be  expressed  as  a closed 
form.  The  forms  given  below  are  readily  evaluated  by  means  of  a 
small  digital  calculator. 


PA  * Q , for  a * 1 


(6.32) 


Q°(b  + 1 - bQ),  for  a ■ 2 


(6.33) 


h Qb  [(b+2)  (b+1)  - 2(b+2)bQ  + (b+l)bQ2],  for  a = 3 (6.34) 


■it  S' 
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Transition  curves  for  six  cases  are  given  in  Figures  6.8 
through  6.10,  which  also  show  the  transition  curve  for  a receiver 
with  a two-channel  threshold  operating  in  Gaussian  noise.  It  is 
seen  that  the  spread  is  primarily  determined  by  the  parameter  a. 
(The  coefficient  of  variation  of  the  signal  power  envelope  is 
1 t -^al ) The  case  a = 1 corresponds  to  a Rayleigh  fading  signal. 


r 


FIGURE  6.9  Transition  Curves  for 


FIGURE  6.10  Transition  Curves  for 
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7.0  SPECIAL  CASE:  P.(t)  * PN(t) 

This  section  Investigates  the  special  case  In  which  the 
noise  processes  have  a common  power  envelope  process:  i.e., 
Pi(t)  * PN(t),  i = 0,  1,  2,  . ..,  n.  This  case  is  relevant 
to  the  case  of  a multi-beam  system  whose  performance  is  limited 
by  wave  noise.  Since  noise  level  is  related  to  average  wind 
speed,  and  since  average  wind  speed  is  correlated  over  a large 
area,  it  seems  reasonable  to  expect  the  power  envelopes  of 
different  channels  to  be  highly  correlated.  A similar  phenom- 
enon might  be  expected  with  a frequency  analyzer  operating 
in  a frequency  band  in  which  the  ambient  noise  is  dominated 
by  shipping.  Since  all  channels  are  viewing  the  same  sources, 
it  is  reasonable  to  expect  a high  correlation  of  noise  power 
envelope  from  bin  to  bin  if  the  power  spectra  of  the  sources 
are  reasonably  smooth.  The  case  to  be  examined  is  a limiting 
case  in  which  the  noise  power  envelopes  are  perfectly  corre- 
lated. 


Additional  assumptions  are  listed  below. 

1.  The  noise  power  envelope  process  is  given 
by  (4.23). 

2.  The  parameters  of  all  channels  are  identical. 

3.  The  signal  is  a stationary  Gaussian  process; 
therefore,  Ps(t)  ■ Ps,  a constant. 

4.  The  weighting  coefficients  c^  for  threshold 
formation  are  those  derived  in  Section  6.1. 
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5.  The  number  of  channels  employed  for  threshold 
formation  is  large. 


6.  All  spectral  functions  have  the  same  rec- 
tangular passband  of  width  W. 


The  initial  objective  of  the  analysis  is  to  evaluate 
the  first  three  moments  of  the  test  statistic  for  this  case. 


For  the  general  case,  the  mean  value  of  the  test  statistic 
is  given  by  (5»3)j  for  the  case  at  hand,  this  may  be  expressed 


PS  ~ ?N  1=0  C± 


(7.1) 


Given  Assumptions  4 and  5,  the  indicated  sum  of  the  weighting 
coefficients  is  y * ^ WT,  where  y is  a false  alarm  parameter, 

W is  input  bandwidth,  and  T is  averaging  time.  Thus 


Ps  '^Pn 


(7.2) 


where  pN  = E[PN(t)] 


For  the  general  case,  the  variance  of  the  test  statistic 
is  given  by  (5.17);  for  the  case  atyhand,  this  may  be  expressed 
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o*  = (WT)"* 


T 

where  In  = T-1  / dy  (1  - | y 1 T” 1 ) r2(y) 

-T 

is  the  variance  of  PN(t),  the  noise  power 
envelope 


r(  ) is  defined  per  (4.30)  or  (4.46) 

n 


Given  Assumptions  4 and  5,  E c?  “ 1:  thus,  after  the  terms 

i-0  1 

are  regrouped,  the  result  is 


o2  - (WT)-1 


(pstp4  +0P2  (1  + ^ b) 


(7.4) 


For  the  general  case,  the  expression  for  the  third  central 
moment  of  the  test  statistic  is  given  in  Table  2,  in  Section  5.0. 
Applying  all  of  the  assumptions  except  5 gives 


continued. . . . 
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* V.  * A A ^ 


n n 


PNeP  + 2pN°P 


P“6p  (ll  °l)  rT 


(7.5) 


where  I*,  * 6T“2  I dx  r(x)  (1  - xT_i)  / dy  r(x-y)r(y) 
1 0 0 


6p  * 23n”2  for  the  power  envelope  process  defined  by 

(4.23) 

* ( 2 ( 1— a ) ] 3 t (1-a3)  for  the  envelope  process 
defined  by  (4.37) 


With  the  application  of  Assumption  5,  the  sum  E is 
approximately  minus  one;  thus  1 0 


- <TW)-J  j «jpNp|  + Ps(pS  ♦ + PNPjJ 

♦ 2 [p|  + Pn(!  + 6p)l 

- 3Y  >[W  ^2^Ps  + Pj|)°lp  + pN6pj  JD 
-(ypnJ  ^Jr  ep  iT[ 


(7.6) 


rMf 


xMH* 
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If  the  square  root  of  the  time-bandwidth  product  is  very 
large,  certain  terms  in  the  expressions  for  the  variance  and 
the  third  central  moment  of  the  test  statistic  may  be  neglected. 
The  ratio  of  the  mean  of  the  test  statistic  to  its  standard 
deviation  is  then 


(7.7) 


where 


r " PS  * PN 


rT  * Y * \]  WT 


B * 1 + Cp  (1  + y2Id) 


CP  = ?N  ‘ aP 


The  coefficient  of  skewness  of  the  test  statistic  is  then 


y3Z  * aZ 


-Y  ^3p  + 2C ^ ID  + y20p  IT]  B-  V2  (7.8) 


The  integrals  Ip  and  I»p  have  been  evaluated  for  the  case  in 
which  the  autocovariance  function  for  the  power  envelope  is 


Kp(u)  - o|  exp  1 1 1 D-1^ 


(7.9) 
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where  D is  the  relaxation  time  of  the  envelope  process.  For 
the  power  envelope  processes  described  in  Section  4.3, 


Kp(u)  = [apr(u) ] 2 


(7.10) 


where  r(  ) is  the  autocovariance  function  for  the  Gaussian 
processes  whose  squares  are  the  components  of  the  envelope 
process.  Using  (7.9)  and  (7.10)  gives 


r(u)  = exp  [-  1 1 1 (2D)-1  ] 


(7.11) 


Evaluating  the  integrals  shown  below  (7-3)  and  (7-5)  with  this 
function  gives 

ID  * (4D  4 T)  [1  + (2D  4 T)  (exp  - T 4 2D  - 1)]  (7-12) 

IT  = 6 (2D  4 T)2  [1  - 4D  4T  + (1  + 4D  4 T)  exp  - T 4 2D]  (7-13) 

The  results  given  above  will  first  be  applied  to  cases  in 
which  the  distribution  of  the  test  statistic  deviates  moderately 
from  normal,  and  (3-11)  can  be  used  to  compute  the  probability 
of  detection.  The  noise  envelope  process  is  assumed  to  be  a 
normalized  chi-square  (NCS)  process.  Figure  7.1  shows  transi- 
tion curves  for  cases  in  which  the  NCS  process  has  four  degrees 
of  freedom,  yielding  a kurtosis  of  4.5  for  the  input  noise  pro- 
cess.* It  is  seen  that  the  spread  of  the  transition  curves 
increases  with  the  ratio  D/T  of  the  envelope  relaxation  time  to 
the  receiver  averaging  time.  The  transition  curve  pertaining 
to  Gaussian  noise  is  also  shown.  Figure  7.2  shows  transition 
curves  for  cases  in  which  the  NCS  process  has  eight  degrees  of 
freedom,  yielding  a kurtosis  of  3.75  for  the  input  noise  process. 


•For  a Gaussian  process,  the  kurtosis  is  3-0. 
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FIGURE  7.1  Transition  Curves  for  NSC 
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Both  figures  show  that  small  values  of  D/T  produce  significant 
spreading  of  the  transition  curves. 


For  larger  values  of  D/T,  the  distribution  of  the  test 
statistic  will  be  considerably  non-normal,  and  (3.11)  cannot 
be  used  to  predict  the  probability  of  detection.  In  order 
to  proceed  further,  an  alternative  to  the  normal  density  is 
required  for  a basis  function  in  (3.4).  A density  function 
for  the  test  statistic  will  be  derived  for  the  limiting  case 
described  in  Section  6.2,  in  which  P,  (t)  = PM;  that  is,  the 
noise  envelope  is  a time-invariant  random  variable. 

In  the  case  discussed  above,  the  multichannel  input 
noise  processes  are  given  by 


Ni(t)  = Gi(t)  (7.14) 


For  the  case  of  noise  alone,  substituting  (7.14)  in  (2.9)  gives 


n T 

Z = -P  E c,  T“  / du  Gl(u)  (7.15) 

i=0  0 1 

This  result  can  be  expressed  as  the  product  of  two  independent 
random  variables: 

B * PW  (7.16) 


where  B = Z for  this  special  case 


n T 

W - - E c,  T-1  / du  G?  (u)  (7.17) 

i-0  1 0 1 
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If  the  bandwidths  of  the. independent  normal  processes  G^(t)  are 
very  large  compared  to  the  reciprocal  of  the  post-rectification 
averaging  time  T,  then  the  variates  represented  by  the  integrals 
in  (7.17)  are  nearly  normal;  furthermore,  a linear  combination 
of  these  variates  is  even  more  nearly  normal.  Thus,  in  most 
cases  of  interest,  the  random  variable  W can  be  considered  as 
normal,  and  the  density  of  B can  be  obtained  if  the  density  of 
P is  specified. 

One  form  for  the  density  function  of  B is 


fR(x)  ■ / dy  fp(y)fw(xy-1 ) y-1 


where  fp(  ) is  the  density  function  for  P;  it  is  zero 
for  negative  values  of  its  argument 


fw(  ) is  the  density  function  for  W. 


An  alternative  form  is 


fB(x)  * / dy  fp(xy-1 )fw(y)  y“l,  x > 0 


(7.18) 


(7-19) 


- / dy  fp(xy-l)fw(y)  y1,  x < 0 (7.20) 


The  use  of  the  results  of  Section  3.0  requires  a basis 
function  for  a standardized  random  variable 


B - m, 


(7.21) 


B m ■ 
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and  the  density  for  this  variable  is 


fu(x)  = VlJ^gX  + %) 


a. 22) 


Substituting  this  result  in  (3.8),  and  that  result  in  turn  in 
(2.12)  gives  the  probability  that  the  test  statistic  exceeds 
the  threshold  value  (zero): 


PE  ’ X0  + (“3  ' m3)  A‘‘  J0  41  X1 
(®BX  + D1b)  xl 


where  1^  - <Jg  / dx  fg  I afix  + mB 

^V°Z 


Let  u = OgX  + mgj  then 


(7.23) 


I = / du  f„(u) 
1 L B 


U-! 


“b 


aB 


(7.24) 


where  L = I — 

B V Og  az 


Substituting  (7.19)  and  (7.20)  in  (7.23)  gives 


1^  » + B^ , L <_  0 


where  A.  - / dy  f^yjy**1  / du  fp 

1 o " 0 r 


B.  - -/  dy  f, 


(-■)(=?) 


(7.25) 


(7.26) 


(7.27) 
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In  the  inner  integrals,  let  v = uy"1;  then 


r-v 


Ai  * l-sf  I ' V*5  ' iv  fp<v) 


(l  - vying') 


(7.28) 


JB  / -» 


Ly-1 
I 
0 


Bi  = l^f)  dy  fw(y)  { dv  fp(v)^l  - vying1  ^ (7. 


29) 


In  the  Integrands,  the  polynomials  can  be  represented  as  sums: 

(7.30) 


^ j A c, 

1 \ By  j-o  w J 


B 


1 • (?)  jo  0 D; 


(7.31) 


where  C, 


w 


00  ^ 

r 00 

o dy  fw(y)yJ  / dv  fp(v)vJ 


(7.32) 


W"J  0 . Ly-1  1 

/ dy  fw(y)yJ  / dv  fp(v)vJ 


(7-33) 


0 


is  the  binomial  coefficient 


In  (7.32),  the  inner  Integral  gives  the  moments  of  the 
random  variable  P;  thus 


ci  • v " dy  fw(y)yJ 


(7.34) 
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If  the  random  variable  W is  assumed  to  be  normal,  then 


fw(y) 


’wV17 


exp  -1/2 


(?J 


(7.35) 


If  (7-35)  is  substituted  in  (7.34),  and  if  the  variable  of 
integration  is  changed  to  x - (y-mw)  * ow  then 


. 00 

C1  " (-Pn)_J  miP  “ t dx  (1  - xy-1)-3  exp  -x2/2  (7-36) 

J J JTtt  y 

since  mg  = Pj^rriy.  If  the  polynomial  in  .the  integrand  is  repre- 
sented as  a sum,  and  if  the  order  of  integration  and  summation 
are  interchanged,  then 


■ <-Pn>-J  mjP  j0  (J)  <-*>' 


(7.37) 


— I dx  x*3  exp  - x2/2 
^ y 

(7.38) 

00 

; — / dx  exp  - x2/ 2 = Q(y) 

^2?  y 

(7.39) 

— exp  -x2/ 2 = n(y) 

^27 

(7.40) 

Yn(Y ) + Q(Y) 

(7.41) 

(2  + y)  n (y) 

(7.42) 
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If  (7.39)  through  (7.42)  are  substituted  as  required  in  (7.37), 
and  if  the  results  are  substituted  as  required  in  (7.30),  the 
end  results  are 


A0  - Q(y)  (7.43) 

A1  * IWb1  n(Y)  (7.44) 

A2  = (pN°WaB1)2  K1  - cp)  Yn(y ) + [l  + C*(l  + y2)]q(y)J 

(7.45) 

A3  * (‘VVi1)  { [6cP  + (2  + Y2)  (l  + W3p*Pi5)]n(Y> 

- ^6 c*  + (3  + Y2  yQ(y)|  (7.46) 

To  evaluate  (7.33),  assume  that  is  a gamma  variable. 
Substituting  (6.26)  in  (7-33)  gives 


■ (-"b)j  i dy  f«<y,yj  ^ Y dv  teR 


(7.47 


Now  let  x • bvp“l;  then 


/ p yJ  0 bL(pNy)“l 

Di  " (isrs)  f dy  fw(y)yJr‘x(b)  / **  xb+J_1  e-* 
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The  inner  integral  is  evaluated  by  means  of  3-381  No.  1.  of 
Reference  13: 


dj*  =4*  L dyf 


w(y)yJ  yjb+j,  bL(pNy)"‘j 


r"1 (b) 


(7.49) 


where  y(  , ) is  the  incomplete  gamma  function.  If  the  random 

variable  W is  assumed  to  be  normal,  and  if  (7.35)  is  substituted 
in  (7-49),  and  if  the  variable  of  integration  is  changed  to  x = 
-y(YOw)“*  the  result  is 

J OB  r j T 

D.  ■ v ^ ~|b ^ / dx  exp  - X (1-x)2  y(b+J,  bKx-1)r-1(b) 

J JET  0 L * J 


where  K ■ -L(pNyow)_l 


(7-50) 


(7-51) 


If  the  quantity  L,  defined  below  (7.24),  is  substituted  in  (7-51)» 
the  result  obtained  after  algebraic  operations  is 


-H-*) 


(7-52) 


Then,  utilizing  (7-2)  and  (7-4)  and  simplifying  the  result  gives 
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K - 1 + 


(r+1) 


\(r+lV 


2 + Cp(l+y2 ) /r  \ 

+ =|(i+vaiD)  \r~ ' 7 


(7.53) 


For  the  small-signal  case  (r  = ps  * < < 1 ) , the  parameter  K 

depends  primarily  on  the  ratio  of  r to  rT. 

If  the  noise  distribution  parameter  b is  an  integer  n, 
then,  by  8.352  No.  1 of  Reference  13, 


Y(b+J,  bKx“ 1 ) 

T(b) 


(b+J-1 ) ! 
(n-1)! 


1 - 


b+J-1  (r\  Kv-1  \m  1 

Z ) exp  -bKx“l 

m-0 


(7.54) 


The  coefficients  required  for  (7.22)  are  derived 

from  the  moments  of  the  standardized  random  variable  U given 
by  (7-21)  using  the  determinants  defined  under  (3*8).  These 
moments  are  derived  from  the  moments  of  the  zero-mean  random 
variable 


V - B - 


PW  - 


“ "PNinW 


1 - PP-1  (1  - 0Y“») 
N 


(7.55) 


where  0 is  a standardized  normal  random  variable.  The  moments 
of  V are  given  by 


f 

{ 


f 

I 


r 


107 


Bolt  Beranek  and  Newman  Inc. 


Report  No.  3656 


E[Vn] 


n /n 


= (-P^)"  Z [ ] m1P(_PN)  1 A I 


1=0  V 1 


J = 0 \ . I mlG 


(7.56) 


If  P Is  a gamma  random  variable, 


m.p  p-i  * ^ 1 

iP  N (b-l)lb1 


(7.57) 


and  the  moments  of  G are 


m21G 


(21)! 
1 ! 21 


m(2i+l)G  = 0 (7.58) 

Values  for  plotting  transition  curves  were  calculated 
using  (7.23)  in  which  the  value  of  the  determinant  A and  values 
of  6i  were  obtained  from  the  determinants  given  below  (3.8). 

The  required  moments  were  calculated  from  (7.56)  using  (7.57) 
and  (7.58).  The  normalized  moment  a ■ + a*  was  obtained 

from  (7.8).  The  integrals  1^  given  by  (7.25)  were  evaluated 
using  (7.43)  through  (7.46)  and  (7.51).  The  latter  were 
executed  by  means  of  a digital  computer. 

The  transition  curves  are  shown  in  Figure  7.3  for  the  cases 
in  which  the  power  envelope  is  a normalized  chi-square  with 
four  degrees  of  freedom.  The  case  T + D * 0 is  that  in  which 
the  envelope  is  a random  variable;  i.e.,  time  invariant. 
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FIGURE  7.3  Transition  Curves  from  Non-normal  Basis. 
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Appendix 

An  Alternative  Approach 

An  alternative  approach  to  that  discussed  in  Section  3.0 
for  predicting  detection  performance  is  based  on  the  assumption 
that  the  averager  outputs,  and  hence  the  test  statistic,  are 
conditionally  normal*  given  the  power  envelope  functions. 

The  probability  of  the  test  statistic  exceeding  zero  is 
given  by  (2.12)  in  terms  of  the  density  of  the  standardized 
test  statistic  Q.  In  analogous  fashion,  the  conditional  prob- 
ability can  be  expressed  as 


P(Z>0|P)  = / dq  fQ(q|P)  (A-l 

"mZC*aZC 

where  P s (P^t)},  i = S,  0,  1,  2,  ...»  n 
mzc  is  E(Z| P) 
a|c  is  Var ( Z | P ) 

If  Z is  conditionally  normal,  then 

P(Z>0|P)  - N(mzc  * ozc)  (A-2 

where  N(  ) is  the  normal  distribution  function  for  a variate 
with  zero  mean  and  unit  variar 


•An  approach  based  on  conditional  norma  *v.y  was  described  in 
Reference  16  and  applied  to  a case  with  a Rayleigh  Fading 
signal  and  normal  noise. 
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An  expression  for  the  test  statistic  Is  given  by  (5.1)* 
If  the  inputs  to  the  multichannel  analog  are  processes  of  the 
type  given  by  (4.5),  the  test  statistic  is  then 


Z = 2T"X  I du  ^Ps(u)PQ(u)  Gs(u)GQ(u) 


n T 

2 c,  T”1  f du  P. (u)G2(u) 
1=S  1 0 1 


(A-3) 


The  mean  value  of  Z conditioned  on  the  power  envelope  functions 
is 


n 


E(Z  P)  = - E c.  A, 
i=S  1 1 


(A -4) 


where  A±  = T"1  / du  P^u)  is  a random  variable  that  is  a 
0 finite  time  average  of  the  1th 

power  envelope. 

The  conditional  variance  of  the  test  statistic  is  the  conditional 
expected  square  less  the  conditional  mean  squared;  the  result  is 


Var  (Z | P)  = 4T"2  / du  / dv  ^Ms(u)Ms(v)M0(u)MQ(v)ps(u-v)p0(u-v) 


n T T 

+ E c2  2T“2  / du  / dv  M, (u)M,  (v)p2(u-v) 
i»S  1 0 0 1 1 


where  p,(  ) is  the  autocovariance  of  G^t), 


(A-5) 
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Convenient  approximations  for  the  terms  of  (A-5)  will  be 
derived.  With  a change  of  variables  x * u+v  and  y = u-v,  the 
1th  integral  of  the  second  term  is 


0 y+2T 

T“2  dy  p£(y)  y dx  P±  [}s(x+y ) ] P.^  teU-y)  ] 

-y 


-T 


T -y+2T 

+ T“2  I dy  p!(y)  / dx  P, [*s(x+y ) ]P, [*s(x-y) ] 

0 1 y 1 1 


(A -6) 


If  the  correlation  time  of  the  process  Gi(t)  is  much  shorter 
than  the  correlation  time  of  the  power  envelope  process,  the 
contributions  of  the  inner  integrals  are  significant  only  in  a 
small  region  around  y » 0;  thus 


T 2T 

I a ij-2  y dy  p.2(y)  / dx  M? (x/2 ) 
1 -T  1 0 


(A-7) 


If  the  correlation  time  of  G^(t)  is  much  shorter  than  T,  and  if 
the  variable  of  integration  of  the  second  integral  is  changed  to 
u * x/2,  then 


I * T"1  W-1  B 
xi  1 l ci 


(A-8) 


where  W"1  * 2 / dy  p2(y) 


B,  * T”"1  / du  P?(u)  is  a random  variable  that  is 
* 0 1 

the  finite  time  average  of  the 
square  of  the  1th  power  envelope. 


114 


I 


Bolt  Beranek  and  Newman  Inc. 


Report  No.  3656 


A similar  procedure  is  employed  for  the  first  term  of  (A-5); 
utilizing  that  result  and  (A-8)  gives 


Var  ( Z | P)  = T_1  (2W£l  Bc  + Z c*  wi*  B±) 

i=S 


(A-9) 


where  W”1  -21  dy  Ps(y)pQ(y) 


Bc  = T~1  f du  Ps(u)PQ(u) 


With  regard  to  the  mean  value,  the  condition  on  the  power 
envelope  process  is  reduced  to  the  condition  on  a set  of  random 
variables  {A^,  i * S,  0,  1,  2,  ....,  n;  and  with  regard  to  the 
variance,  the  set  of  random  variables  is  {Bj } j ® C,  S,  0,  1,  2, 

3,  ....,  n.  The  combined  set  is  a random  vector 

C - {{A^,  {Bj}} 

For  the  general  case,  it  is  seen  that  the  number  of  random  variables 
can  be  large;  the  exact  number  depends  on  the  number  of  non-zero 
weighting  coefficients  c^.  Furthermore , and  B^  are  always 
dependent,  as  are  Ag,  AQ,  Bc,  Bg,  and  BQ.  And  in  some  cases,  the 
sets  {A^}  and  {B^>  could  be  dependent. 


The  probability  that  the  test  statistic  exceeds  zero  is  the 
expected  value  of  the  conditional  probability 


/ ....  f dcfc (c)  N 
0 0 


E[Z|c] 


(A-10) 
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where  f ( ) is  the  Joint  probability  density  function  for  the 
c 


component  of  C,  and  where  E [ Z | ] , Var  [Z|_]  are 

aiven  bv  (A-4)  and  (A-9)  resDectivel.v. 
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